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ABSTRACT 

Circumstellar disks are thought to experience a rapid "transition" phase in their evolution that can 
have a considerable impact on the formation and early development of planetary systems. We present 
new and archival high angular resolution (O'.'S « 40-75 AU) Submillimeter Array (SMA) observations 
of the 880 /zm (340 GHz) dust continuum emission from 12 such transition disks in nearby star-forming 
regions. In each case, we directly resolve a dust-depleted disk cavity around the central star. Using 
two-dimensional Monte Carlo radiative transfer calculations, we interpret these dust disk structures 
in a homogeneous, parametric model framework by reproducing their SMA continuum visibilities and 
spectral energy distributions. The cavities in these disks are large (i?cav = 15-73 AU) and substantially 
depleted of small (~/zm-sized) dust grains, although their mass contents are still uncertain. The 
structures of the remnant material at larger radii are comparable to normal disks. We demonstrate that 
these large cavities are relatively common among the millimeter-bright disk population, comprising 
at least 1 in 5 (20%) of the disks in the bright half (and >26% of the upper quartile) of the millimeter 
luminosity (disk mass) distribution. Utilizing these results, we assess some of the physical mechanisms 
proposed to account for transition disk structures. As has been shown before, photoevaporation models 
do not produce the large cavity sizes, accretion rates, and disk masses representative of this sample. It 
would be difficult to achieve a sufficient decrease of the dust optical depths in these cavities by particle 
growth alone: substantial growth (to meter sizes or beyond) must occur in large (tens of AU) regions 
of low turbulence without also producing an abundance of small particles. Given those challenges, we 
suggest instead that the observations are most commensurate with dynamical clearing due to tidal 
interactions with low-mass companions - young brown dwarfs or giant planets on long-period orbits. 
Subject headings: circumstellar matter — protoplanetary disks — planet-disk interactions — planets 
and satellites: formation — submillimeter: planetary systems 



1. INTRODUCTION 

Planets are made from gas and dust in disks that or- 
bit young stars. The physical characteristics and evo- 
lution of that disk material have a profound impact on 
the outcome and efficiency of the planet formation pro- 
cess. While a massive theoretical effort has been made 
to elucidate the physical mechanisms that contribute to 
the assembly of a planetary system, any successful planet 
formation model needs to be informed by observations. 
An important goal is to generate models that incorpo- 
rate empirical insights on the properties of circumstellar 
disks and are able to reproduce ke y demographic features 
of the exoplanet population ('e.g..llda fc Linll2004 120051 : 
lAlibert et all 120051: iMordasini et al.ll2009[) . One essen- 
tial aspect of meeting that goal is the ability to make 
measurements related to the initial conditions imposed 
by the disk material. However, those constraints on ini- 
tial conditions must be coupled with a characterization 
of disk evolution, which is ultimately responsible for the 
dynamical and physical properties of any resulting plan- 
etary system. 
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For most of its lifetime, the evolution of a protoplan- 
etary disk is regulated by the interaction of gravita- 
tional and viscous stresses (e.g.. lPringlelll98ll ). Thought 
to be produced by ma gnetohydrodynamic turbulence 
(|Balbus fc Hawlevlll99l . the viscosity of the disk ma- 
terial acts to channel mass toward the star by trans- 
porting angular momentum o ut to larger dis k radii 
(|Lvnden-Bell fc Pri^j^ HqtI iHartmann eFall [l998.) . 
This viscous evolution process decreases the densities and 
mass flow rates in the region of the disk where plan- 
ets are made, and therefore sets a bound on the period 
where formation conditions are most favorable. More- 
over, the viscous properties of the disk control the dy- 
namic s of solid bodies, from the radial drift of dust par- 
ticles (lWeidenschilling||1977aD to the o r bital migration of 
protop lanets ()Lin fc Papaloizoulll986l ). IHartmann et al.l 
(|1998f) introduced some indirect methods for placing ob- 
servational constraints on viscous evolution, based on 
tracking the decline of a ccretion rates with age (e.g., 
iSicilia-Aguilar et "al1l2010t ) and characterizing disk sizes 
and surface density profiles with resolved millimeter ob- 
servations (iKitamura et al.l 120021: [Andrews et"all 120091 
I2010btllsen"a et al.ll2009ll2010aD . More direct approaches 
have focused on characterizing the anomalous viscosities, 
including attempts to measure the disk magnetic field 
structure d Hughcs et al.l 120090 ) and turbulent velocities 
(jHughes et al..,2011i) . 

Although viscous evolution controls the long-term be- 
havior of disk material, there is strong observational evi- 
dence that some (if not all) disks undergo a faster evolu- 
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tionary phase that should have important consequences 
for the assembly and development of planetary systems. 
A small subsample of young disks exhibits a distinctive 
dip in their infrared spectral energy distribution (SED), 
suggesting that warm dust near the central star is sub- 
stantially depleted compared to larger ra dii in the disk 
(jStrom et al.lll989t iSkrutskie et al.l ll990[) . Since these 
so-called "transition" disks represent only a few per- 
cent of the total population in nearby star-forming re - 
gions with ages up to a few Myr (|Muzerolle et al.ll201Cll ). 
the inside-out dispersal of disk material they appear to 
trace must occur quite rapidly (a few x 0.1 Myr, as- 
suming all disks go through this phase at some point). 
A variety of dissipation mechanisms have been pro- 
posed to explain observations of t ransition disks, in- 
cluding extensive part i cle gr owth ([Tanaka et al.l I2Q05I : 
IDullemond fc DominikI I2005D. photoevapor ative mass - 
loss in winds HClarke et al.ll200 It [Alexander et"aIll20Q6aD. 
and tidal interactions with co mpanions (e.g.. iRice et al.l 
I2003at [Ireland fc KrausI 120081) . As a tantalizing exam- 
ple of the latter possibility, the transition disk struc- 
tures might have been modified by dynamical interac- 
tions with very young giant planets (e.g., Bryden et ahl 
ll999HPapaloizou et al.lf2007i ). If that is the case, the size 
and material content of the dust-depleted region can in 
principle be used to constrain the o rbit and mass of the 
unseen planetary compan ion (e.g., iLubow fc D'Angelol 
[200l IVarniere et al.l[2006h . 

Regardless of which physical mechanism is responsi- 
ble, this rapid evolutionary phase is potentially extraor- 
dinarily important in the context of planet formation. 
Since the short timescale for the dispersal or transfor- 
mation of disk material can be quite influential on the 
accretion and migration history of young planets, this 
transitional period may play a dominant role in setting 
the final masses and orbital configurations of planetary 
systems. In that sense, the transition disks are crucial 
benchmarks for planet formation models that rely on ob- 
servational comparisons of disk properties with the de- 
mographic characteristics of the much older population of 
exoplanets. There is an extensive literature on transition 
disk observations, ranging from their identification and 
initial probes of e i npirical trends (e .g..lNaiita et al.ll2007l : 
iCieza et al.|[2008l: iKim et al.|[200a iMerm et al.ll2010D to 
detailed analyses of their b r oadband SEDs and infrared 
spectra (Calvet et al.' '200^, '2005|; ID'Alessio etHI 12001 
lEspaillat ct al. 2007, 2010, 201l5. To facilitate the in- 
clusion of this phase of disk evolution into the next gen- 
eration of planet formation models, such studies must 
be coupled with firm observational constraints on the 
sizes and material contents of the dust-depleted zones in 
transition disks, as well as improved characterizations of 
their remnant mass reservoirs at larger radii. Those goals 
are best achieved by resolving their optically thin dust 
emission at millimeter wavelengths (Piitu et al.l 120061 : 
iBrown et al.][2008l [20091: iHughes et al]|2009aD . 

In this paper, we present the first large sample of high 
angular resolution millimeter continuum observations of 
protoplanetary transition disks and analyze them in a ho- 
mogeneous framework by constructing parametric mod- 
els of their dust structures. These interferometric data 
and their calibration are described in §2. Detailed ex- 
planations of the adopted modeling formalism, radiative 
transfer calculations, and uncertainties in the modeling 



process are provided in §3. A preliminary characteriza- 
tion of the transition disk structures is included in §4. 
The data and models are synthesized and discussed in 
the context of disk evolution mechanisms in §5, with the 
key results summarized in §6. A brief commentary on 
each individual disk is included in the Appendix. 

2. OBSERVATIONS AND DATA REDUCTION 

We have constructed a sample of 12 nearby transition 
disks with high angular resolution 880 fim continuum 
data using new and archival ob servations with the SMA 
interferometer ()Ho et al.ll2004D . These targets were se- 
lected either based on the presence of a distinctive dip in 
their Spitzer IniraRed Spectrograph (IRS) spectrum (see 
the Appendix for references) and bright 880 /ini emis- 
sion (Fi, > 100 mjy) or the previous confirmation of 
a large dust-depleted cavity via direct imaging. Those 
criteria introduce some important biases that are dis- 
cussed in more detail in §4. A journal of the SMA ob- 
servations is provided in Table [1] Since some of the 
data in this sample w e re prev i ously described els e where 
JAndrews et al.' '^OOg', '2010 bl: IBrown et al.l [200l l200l 
iSughcs ct al. 2009a; IscUa et al.ll2010bO . details of the 
adopted instrument configurations and calibration pro- 
cess will not be repeated here. New observations of DM 
Tau, LkCa 15, UX Tau, and RX J1615.3-3255 were ob- 
tained in 2010 February with the "very extended" (V) 
configuration of the SMA, with baseline lengths ranging 
from 68 to 509 m. Those data were supplemented with 
shorter baseline observations in the compact (C: baseline 
lengths of 6-70 m) and extended (E: baseline lengths of 
28-226 m, excluding RX J1615-3255) configurations over 
the past year (since 2009 October), as were additional ob- 
servations of GM Aur (C, E), LkHa 330 (C), and MWC 
758 (C). 

For these recent observations (since fall of 2009), we 
used an expanded bandwidth capability (with the ex- 
ception of the C observations of DM Tau) that sepa- 
rates each sideband into two intermediate frequency (IF) 
bands spanning 3-5 GHz and 5-7 GHz from the local os- 
cillator (LO) frequency (340.755GHz, or 880 ^im). Each 
IF band is composed of 24 partially overlapping 104 MHz 
chunks, with the central chunk in the lower IF band split 
into 512 channels to sample the CO J=3— 2 transition 
(345.796GHz) at a resolution of 0.18km s'^ All other 
chunks were divided into 32 coarse channels to measure 
the continuum. Using a cycle time of ^10 minutes, the 
observing sequence interleaved disk targets with nearby 
quasars: J0510-f 180, 3C 111, and 3C 120 for the north- 
ern disks and J1626-298, J1625-254, and J1604-446 for 
RX J1615-3255. When the targets were at low eleva- 
tions (<20°), planets (Uranus, Mars), satellites (Titan, 
Callisto), and bright quasars (3C 273 and 3C 454.3) were 
observed as bandpass and absolute fiux calibrators, de- 
pending on their availability and the array configuration. 
Weather conditions were excellent, with 225 GHz atmo- 
spheric opacities below 0.1 (2.0 mm of precipitable water 
vapor) in all cases and substantially lower (down to 0.03, 
or 0.6 mm of water vapor) for some observations. 

The data were calibrated using the MIR software pack- 
age, treating the individual IF bands separately to in- 
dependently assess the accuracy of the process. The 
bandpass response was corrected based on observations 
of a bright quasar, and broadband continuum channels 
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Table 1 

SMA Observine Journal 



Name 


a [J2000] 


(5 [J2000] 


Region 


Array 


UT Date 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


MWC 758 


05 30 27.53 


+25 19 56.9 


Tau? 


E 


2008 Jan 28 


1 










C 


2010 Nov 3 


2 


SAO 206462 


15 15 48.43 


-37 09 16.3 


SCO-OB2? 


V 


2007 May 27 


3 










V 


2007 Jun 8 


3 


LkHa 330 


03 45 48.29 


+32 24 11.9 


Per 


V 


2006 Nov 19 


4 










c 


2010 Nov 3 


2 


SR 21 


16 27 10.28 


-24 19 12.8 


Oph 


V 


2007 Jun 10 


3 










s 


2008 Aug 29 


5 


UX Tau A 


04 30 04.00 


+18 13 49.3 


Tau 


c 


2009 Oct 23 


2 










E 


2010 Jan 9 


2 










V 


2010 Feb 5 


2 


SR 24 S 


16 26 58.51 


-24 45 37.0 


Oph 


V 


2009 Mar 25 


6 










C 


2009 May 4 


6 


DoAr 44 


16 31 33.46 


-24 27 37.4 


Oph 


c 


2008 May 13 


5 










V 


2008 Apr 3 


5 


LkCa 15 


04 39 17.79 


+22 21 03.2 


Tau 


c 


2009 Oct 23 


2 










E 


2010 Jan 13 


2 










V 


2010 Feb 6 


2 


RX J1615-3255 


16 15 20.20 


-32 55 05.1 


Lup 


V 


2010 Feb 16 


2 










C 


2010 May 1 


2 


GM Aur 


04 55 10.98 


+30 21 59.3 


Tau 


V 


2005 Nov 5 


7 










c 


2009 Nov 6 


8 










E 


2010 Jan 16 


8 


DM Tau 


04 33 48.74 


+18 10 09.7 


Tau 


V 


2010 Feb 17 


2 










E 


2010 Sep 21 


2 










c 


2010 Oct 7 


2 


WSB 60 


16 28 16.51 


-24 36 58.3 


Oph 


c 


2008 May 13 


5 










V 


2008 Apr 3 


5 



Note. — Col. (1): Name of host star (in order of the stellar spectral type; see 
§3.3 and Table E]). Cols. (2) & (3): Disk center coordinates (see §2). Col. (4): Star- 
forming region. Col. (5): SMA configuration; V = very extended (68-509 m baselines), 
E = extended (28-226 m), C = compact (16-70 m), and S = subcompact (6-70 m). 
Col. (6): UT date of observation. Col. (7): Origina l refer ences for SMA data: [1] - 
llsella e t al' f 2010b'), [2] - this p aper. [3] -IBrown et all (120091 ') . [4] -IBrown ct al. (20081), 
[5] -lAndrews et al. (2009), [6] - I Andrews et all (|2010bl ). [7] - iHuehes etal.. (,2009a ). [8] 
- Hughes et al. [in preparation). 



in each sideband and IF band were generated by av- 
eraging the central 82 MHz of each chunk (except the 
one containing the CO J=3— 2 transition). The visibil- 
ity amplitudes were set based on observations of plan- 
ets or satellites, with a systematic uncertainty of ^10%. 
The antenna-based complex gain response of the sys- 
tem as a function of time was determined with refer- 
ence to the nearest local quasar projected on the sky, 
typically 3C 111, J0510-hl80, or J1626-298. The other 
quasars observed in each track were used to check the 
quality of phase transfer. Based on that information, we 
estimate that the millimeter "seeing" generated by at- 
mospheric phase noise and any baseline errors is small, 
<0'.'l. The calibrated visibilities from different IF bands, 
sidebands, and array configurations were found to be in 
excellent agreement and combined. Because of the time 
elapsed (4-5 yr) between observations of GM Aur and 
LkHa 330, the new da ta were corrected for proper mo- 
tion (Ducourant et al.][2b05) before combining the visi- 
bilities. The CO line data will be presented elsewhere. 

The combined continuum visibilities were Fourier in- 
verted, deconvolved with the CLEAN algorithm, and re- 
stored with a synthesized beam using the MIRIAD soft- 
ware package. Maps of the 880 /im emission toward each 
target are exhibited together in Figure [1] These images 
were synthesized to accentuate structures on small an- 
gular scales. Given the variety in the available Fourier 



Table 2 

Aperture Synthesis Image Statistics 



Name 


i^tot 


^rcc 


-^pcak 


(T 


db 


PAb 




[Jy] 


[Jy] 


[mjy] 


[mJy] 


["] 


[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



MWC 758 


0.18 


0.16 


59 


1.1 


0.82 X 0.55 


93 


SAO 206462 


0.62 


0.34 


54 


3.5 


0.50 X 0.24 


9 


LkHa 330 


0.21 


0.21 


28 


2.1 


0.31 X 0.27 


84 


SR 21 


0.40 


0.26 


59 


2.6 


0.42 X 0.30 


18 


UX Tau A 


0.15 


0.13 


38 


1.5 


0.31 X 0.28 


52 


SR 24 S 


0.55 


0.25 


68 


1.8 


0.37 X 0.26 


14 


DoAr 44 


0.21 


0.21 


62 


2.7 


0.60 X 0.38 


23 


LkCa 15 


0.41 


0.39 


40 


0.9 


0.41 X 0.32 


64 


RX J1615-3255 


0.43 


0.30 


76 


2.8 


0.52 X 0.26 


16 


GM Aur 


0.64 


0.56 


65 


3.5 


0.32 X 0.27 


38 


DM Tau 


0.21 


0.15 


47 


1.3 


0.47 X 0.31 


33 


{inset) 




0.02 


13 


1.3 


0.37x0.21 


31 


WSB 60 


0.25 


0.25 


113 


3.5 


0.59 X 0.39 


25 


{inset) 




0.20 


67 


3.8 


0.44x0.23 


23 



Note. — Col. (1): Name of host star. Col (2): Total 
integrated flux density. Col. (3): Recovered integrated flux 
density in the synthesized maps in Figure [1] Col. (4): Peak 
flux density (in mJy beam~^). Col. (5): RMS (If) noise level 
in a 10"-wide square centered on the stellar position. Col. (6): 
Synthesized beam dimensions. Col. (7): Synthesized beam 
major axis orientation (east of north). 
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Figure 1. SMA aperture synthesis maps of the 880 continuum omission from this sample of transition disks. Each panel is 2'.'7 on 
a side (offsets are referenced to the disk centers listed in Table [T] see §2) and contains a 50 AU projected scale bar in the lower right for 
reference. Contours are drawn at 3 cr intervals, and the synthesized beam dimensions are marked in the lower left corner (RMS noise levels 
and beam dimensions are provided in Table [2]l. The inset images for the DM Tau and WSB 60 disks were synthesized with higher angular 
resolution, and are shown to scale. 



coverage across the sample, a range of visibility weight- 
ing schemes were adopted to make the maps: the typ- 
ical Briggs robust parameter values were 0-1, although 
in some cases uniform (DM Tau, WSB 60) or natural 
(SAO 206462) weighting schemes were used. The syn- 
thesized beam sizes are marked in the lower left of each 
panel and a projected 50 AU scale bar is shown on the 
lower right. The disks located at southern declinations 
have elongated synthesized beams because their low ob- 
serving elevations prevent sufficient long-baseline Fourier 
sampling in one dimension. The same is true for some 
northern sources that were scheduled for only a partial 
observing track (e.g., DM Tau). The image statistics are 
summarized in Table [21 including the total integrated 
flux densities, the flux densities recovered in the syn- 
thesized images in Figure [T] (note that 1 — Frcc/^tot in- 
dicates how much emission is spatially filtered for the 
adopted visibility weighting scheme), peak fluxes, noise 
levels, and beam dimensions. Each target exhibits the 
continuum morphology expected for a resolved emission 
"ring" , with substantially reduced intensities near the 
stellar position. For intermediate viewing inclinations 
(or elongated synthesized beams), these limb-brightened 



rings display a distinctive double-peaked structure on ei- 
ther side of the stellar position. While that peak-to-peak 
separation offers a preliminary indication of the inner 
ring radius, more detailed modeling is required to accu- 
rately infer the disk structures. 

With these unique 880 /im emission morphologies, 
defining a disk center that can be used as a reference 
in such modeling can be problematic. Emission from the 
stellar photosphere is well below our detection limits at 
these long wavelengths, so there is no point-like reference 
in the data. We made an initial estimate of the disk cen- 
ters based on the optical/near- infrared stellar positions, 
corre cted for any known proper rn otion (IPerrvman et aD 
1997HH0g et al.lll998HCutri et al.ll200llDlicourant et all 



20051) . Aside from the multiple systems UX Tau and 
SR 24 (for which proper motion estimates are particu- 
larly uncertain), those positions were near the apparent 
emission centers in the synthesized 880 /im images and 
well within the absolute SMA astrometric uncertainties 
(~0'.'l). Those positions were then refined using a grid 
search that minimized the imaginary component of the 
visibilities (i.e., the position for which the RMS of the de- 
viations from zero imaginary flux was smallest). In prin- 
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ciple, the imaginary component of a symmetric source 
should be zero. Since our models will assume axisym- 
metric structures, this assumption for centering seems 
appropriate. An a posteriori comparison of the data and 
our modeling results (see §4) confirms that there is no 
statistically significant evidence for asymmetric emission 
structures on the angular scales probed by the SMA. The 
adopted reference centers are listed in Table [T] 

3. MODELS OF TRANSITION DISK STRUCTURES 

The physical structures of transition disks are quite 
complex, and the few studies aimed at their characteriza- 
tion using spatiall y resolved data have adopted very dif- 
ferent approaches (iPietu et al. 2006; Hughes et al. 2007, 



2009 



2010 



Brown et al.' ^OOS', 2009t lAndrews et all l2009i 
B; Usella et al, 2010a,^). In this study, our goal is 



to analyze this large sample of transition disks in a homo 
geneous framework, emphasizing the new constraints on 
basic structural parameters that are responsible for the 
key features of the high angular resolution SMA data. 
These initial results should help guide more in-depth, 
individualized studies in the future. The modeling per- 
formed here is based on the previo us work descr i bed b y 
lAndrews et all (pOOl [2010al lH) and lHughes eTall poTol ). 
but includes some modifications meant to imitate (para- 
metr ic ally) the detailed treat ments of Calvct et al. (200 
2005 1) . ImAlessio et al.l (pOOSh. andlEs paillat et all (|200' 



20101 ). In the following sections, we describe a parametric 
model of the disk structure (§3.1), our assumptions about 
the dust populations (§3.2) and stellar properties (§3.3), 
and the radiative transfer calculations and our approach 
for comparing the models with observations (§3.4). Fi- 
nally, we discuss the observable effects of the model pa- 
rameters and important uncertainties (§3.5), and provide 
a concise summary of some key issues related to the mod- 
eling (§3.6) for those less interested in the details. 

3.1. Parametric Density Model 

The model structure includes two distinct radial zones: 
a dust-depleted "cavity" {R < Rcav) and an "outer 
disk" reservoir (i? > Rcav)- We define a global surface 
density profile, appropriate for a simple accretion disk 
with a power-law radial dis tribution of time-independent 
viscosities (ly oc R^; see iLvnden-Bell fc Pringld HQTI 
iHartmann et al.lll998D . 



- ( 



exp 



2-7' 



(1) 



where Rc is a characteristic scaling radius. This global 
profile is applied to the material in the outer disk; 



S(i? > Rc 



Tig . The densities in the cavity are scaled 



down by a constant depletion factor, Scav: such that 
S(i? < i?cav) = Scav^g- To facilitate a more homoge- 
neous comparison of Scav values for a sample with a range 
of cavity sizes (Rcav), we decomposed the cavity into a 
depleted "inner disk" of fixed extent (out to i? = 10 AU) 
and an empty "gap" between 10 AU and i?cav (the arbi- 
trary choice of the inner disk extent is discussed in §3.5). 
Unless noted otherwise, the inner disk is truncated at the 
dust sublimation radius , Rguh ~ 0.07(L^/Lf7))^/^ AU (as- 
suming Tsub = 1500 K: iD^lemond et al.|[200l . Figure 
[2^ shows a generic surface density model for reference. 



The vertical distribution of dust has been modified 
from our previous models, in an effort to more realis- 
tically treat the irradiation heating of the disk atmo- 
sphere. We have introduced a vertical gradient in the 
dust size distribution, such that larger grains are con- 
centrated toward the midplane and smaller grains are 
distributed to larger vertical heights. In practice, this pa- 
rameterization crudely mimics dust sedimentation (e.g., 
iDullemond fc Domimk|[2004bl:lD'Alessio et al.|[2006l) . At 
each radius, a population of "small" grains represents 
a small fraction of the total column density, (1 — /)S, 
and is distributed vertically like a Gaussian with scale 
height h = hc{R/ Rc)^ ■ Conversely, most of the total 
column (/E) is composed of a "large" grain population, 
condensed toward the midplane with a scale height xh 
(where x ^ !)■ I'^ a spherical coordinate system with az- 
imuthal and mirror symmetry, the two-dimensional den- 
sity structures for each dust population are 



Pi 



2TTRh 



^Rxh 



exp 



exp 



7r/2 - e 



7r/2-e 
Xh 



(2) 



(3) 



where the subscripts 's' and 'V denote the small and large 
dust populations, respectively, Q is the (vertical) latitude 
coordinate measured from the pole {Q = 0) to the equa- 
tor (the midplane, 9 = tt/2), and the value of E depends 
on the radial location as outlined above. Note that in 
this coordinate system, the scale height h is an angle: 
the scale height in physical (distance) units is « hR. 
A generic vertical density profile with this parametric 
form is shown in Figure |2]d. 

The distinctive infrared spectra of transition disks re- 
quire small, but important, local modifications to the 
vertical structure n e ar the disk edges a t Rsuh and Rcav 
(iNat ta et al.1 120011: iCalvet et all 120021: lEspaillat et"all 
2010). At those locations, a large surface area of ma- 
terial is exposed directly (or with relatively low atten- 
uation; see §3.5) to stellar irradiation. That exposure 
substantially heats a thin layer of material and "puffs" 
up th e local dust distributi on (jPuUemond &: DominikI 
[2004at iD^essio et al.l l2005l) . To mimic the emission 
signatures of these structures — referred to as the in- 
ner "rim" (at -Rsub) and cavity "wall" (at i?cav) — we 
permit the local value of h to be scaled up parametri- 
cally (to /iriin and /iwaii, respectively). Those perturba- 
tions are joined to the global scale height distribution 
exponentially over a small radial width (fixed here to 
AR ^ 0.1 AU). 

This structure model has 11 parameters: 5 describe 
the surface density profile, {Ec, 7, Rc, Scav, Rcav}, and 
6 others characterize the vertical distribution of dust, 
{he, ip, X, /, ft-rimj ^waii}- Some of thcsc parameters 
can not be constrained with the available data, and were 
therefore fixed to representative values. The surface den- 
sity gradient was kept at 7 = 1, in line with the results 
for "normal" disks with con tinuous dust distributions 
(jAndrews et al.l 120091 l2010b[ ). We set the settling pa- 
rameters to X = 0.2 (the large grains are distributed to 
20% of the scale height, h) and / = 0.85 (85% of the 
total column is composed of the large grain population). 
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Figure 2. (a) A generic surface density model. The outer disk {R > -Rcav) follows the global Eg profile, overlaid as a dashed curve. The 
surface densities in the cavity (R < Rcav) are scaled down by a factor (5cav. The dust-depleted inner disk is truncated at the sublimation 
radius (Rsuh) and extends outward to an arbitrary fixed radius (10 AU). The gap between 10 AU and i?cav is assumed to be empty, (b) A 
generic vertical density profile. The blue curve tracks the "large" dust population, concentrated toward the midplane with a scale height x'l 
(in the models used here, x is fixed to 0.2). The orange curve marks the "small" dust population that dominates in the disk atmosphere, 
and the dashed curve traces the composite profile. The fraction of the column density contributed by the large dust population is denoted 
/ (in the models used here, / is fixed to 0.85). 



The remaining parameters are allowed to freely vary, al- 
though substantial degeneracies remain (see §3.5). 

3.2. Dust Populations 

A parametric density structure must be populated with 
dust grains of a given size distribution and composi- 
tion. For the sake of homogeneity, we assume dust grains 
with miner al abundances as in the inter stellar medium 
(ISM) (see iWeingartner fc Drain j |20Q1[ ). despite some 
evidence in individual cases for different compositions. 
Each dust population has a power-law distribution of 
sizes (s), n(s) oc s~p, from Sj^in = 0.005 /im to a spec- 
ified Smax- The dust properties in the outer disk are 
fixed: "small" grains have s^ax = 1 "large" grains 
have Smax = 1mm, and p = 3.5. The opacity spectrum 
for a given size distribution was derived from Mie cal- 
culations. Since most of the mass is in "large" grains, 
the selection of Smax ^ A = 1 mm is significant because 
it tends to maximize the millimeter-wave dust opacity. 
Therefore, the densities inferred from the optically thin 
emission we have observed with the SMA can be con- 
sidered lower bounds if much larger particles are present 
(|D'Alessio et all 120061: IDraind 120061) . The 880 /im dust 
opacity for the large grain population is 3.6 cm^ g~^- The 
dust size distribution in the inner disk, inner rim, and 
cavity wall are assumed to be identical, but are char- 
acterized by free parameters for their size distribution, 
{smax, p} (see §3.4 regarding the selection of these pa- 
rameter values). Some representative opacity spectra are 
shown in Figure [31 

3.3. Stellar Properties 

Passive irradiation is the only heating mechanism con- 
sidered here, and therefore the properties of the central 
star are pivotal in setting the thermal structure of a 
model. Careful selection of an input stellar spectrum is 
particularly important for modeling transition disks, as 
the infrared disk-to-star contrast is often low. Therefore, 
the details of the stellar spectrum can have a pronounced 
impact on how the small amount of dust inside the cavity 
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Figure 3. The opacity spectra for representative dust popula- 
tions. The opacities for the small (orange) and large (blue) dust 
populations in the outer disk are fixed; the latter is similar to the 
canonical Bcckwitli ct al. (1990) values (gray dashed line). As il- 
lustrative examples, the opacities for two different dust populations 
we have used inside the disk cavities are also shown (green): {smax, 
p} = {10 Aim, 3.5} (solid) and {1 fizn, 2.5}(dotted). 



is interpreted. To derive stellar input spectra, we first as- 
signed effective temperature (Toff) values based on spec- 
tral classificat ions in the literature, using t he co nversions 
advoc ated by iKenvon fc HartmannI (|1995[ ) and iLuhmarJ 
(|1999f l for types earlier or later than MO, respectively. 
Luminosities (i*) were determined by matchin g red- 
dened , scaled spectral synthesis models (Leieune et al.l 
119971 for a fixed Toff) to the optical/near- infrared pho- 
tometry (see the Appendix). We utilized the extinction 
curve derived bv lMcClure (I2009D . which at these wave- 
lengths is identical to the lMathisI (|1990[ ) law for a large 
total-to-selective extinction value, Ry = 5. For cases 
without dynamical mass constraints, this process was it- 
erat ed for diffe r ent st ellar masses (M*) estimated from 
the iSiess et all (|2000} ) pre-main sequence models. The 
stellar parameters used to generate input spectra are 
compiled in Table |3l 
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Table 3 

Stellar and Accretion Properties 



Name 


SpT 




Ay 


d 


L. 


R* 


M, 


M 




logLx 


Ref 






[K] 




[pc] 


[L©] 






[Mq yr- 


] 


r — ll 
[erg s ^] 




(1) 


(2) 




f A\ 

W 


(5) 




[') 




/'n^ 




(10) 


(.11) 


MWC 758 


A8 


7580 


0.0 


200 


15 


2.25 


1.8 


1 X 10" 


8 




1,2,3, ■•• 


SAO 206462 


F4 


6590 


0.3 


140 


7.8 


2.15 


1.6 


6 X 10- 


9 


29.7 


4, 5, 6, 6 


LkHa 330 


G3 


5830 


3.1 


250 


15 


3.75 


2.2 


2 X 10- 


9 




7,5,8, ■■■ 


SR 21 


G3 


5830 


6.3 


125 


10 


3.15 


2.0 


< 1 X 10 


-9 


30.0 


9, 5, 10, 11 


UX Tau A 


G8 


5520 


1.9 


140 


3.5 


2.05 


1.5 


1 X 10" 


8 


30.3 


12, 5, 12, 13 


SR 24 S 


K2 


4990 


7.0 


125 


4.0 


2.70 


2.0 


1 X 10" 


8 


30.0 


14, 5, 10, 15 


DoAr 44 


K3 


4730 


2.2 


125 


1.4 


1.75 


1.3 


9 X 10" 


9 


29.9 


12, 5, 12, 16 


LkCa 15 


K3 


4730 


1.7 


140 


1.2 


1.65 


1.01 


2 X 10- 


9 


< 29.6 


12, 17, 18, 19 


RX J1615-3255 K5 


4350 


0.4 


185 


1.3 


2.00 


1.1 


4 X 10- 


10 


30.4 


20, 5, 5, 21 


GM Aur 


K5 


4278 


1.1 


140 


1.0 


1.85 


0.84 


1 X 10- 


8 


< 29.7 


22, 23, 18, 19 


DM Tau 


Ml 


3705 


0.6 


140 


0.3 


1.25 


0.53 


6 X 10" 


9 


< 29.7 


24, 17, 18, 25 


WSB 60 


M4 


3270 


3.7 


125 


0.2 


1.40 


0.25 


2 X 10" 


9 




26,5,10, ■■■ 


Note. — 


Col. (1): 


Name of host star. Col. (2): 


Spectral type. 


Col. (3); Visual extinction. 



Col. (4): Effective temperature. Col. (5): Luminosity. Col. (6): Estimated distance to star-forming 
region. Col. (7): Stellar radius. Col. (8): Stellar mass, determined either directly from CO spec- 
tral images or indirectly from pre-main sequence evolution models (see §3.2). Col. (9): Accretion 
rate. Col. (10): X-ray luminosity or 3 cr upper limit. Note that the value for UX Tau A may in- 
clude contributions from the additional stars in the multiple system. To our knowledge, no X-ray 
measurements of WSB 60, LkHa 330, or MWC 758 have been published to date. Col. (11): Liter- 
ature references for spectral_classifications, M«, M, and Lx- [1] - IBeskrovnava et al.l 11999"), [2] - 
[Chapillon ct al. (2008), [3] -JEisncr ct al. (200Sl), [4] - IDunkin et all ll 19971). [5] - this paper (sec §3.3), 
6] - iGradv et all (120091) [7] HCohcn & Kuhi (1979), [8] -'Salvk ct al. (2009), [9] - Prato ct al. (2003|V 
lOj -INatta et all (|2006 | ). [11] - Grosso ct al. (2000), [12] - EspaiUat ct al. (2010), [13] - Konis ct al] 
Il200ll). [14] -ILuhman fc Ri ekc (1999). [15] - Pillittcri ct al. (2010), [16] 
IPietu et al.( ( |2007l ). [18"i~ringlcbv ct al. (2009), [19] - Ncuhauscr ct al. 

1997) , [22] -jW hitc & Ghcz (2001), [2 3] - luutrev et ai.l 
Damiani et al (1995), [26] - Wilking et all (|2005[ ). 



(fT999t). [21] - Krauttcr ct al, 
IKenvon fc Hartmann (1995 ), [25] 



- Montmcrlc et al. ( 198 3)), [17] 
( 1 995,), [20] - Wi c hmann et all 
- IDutrev et al.l (jT998h . 



[24] 



3.4. Radiative Transfer and Modeling Procedure 

For a given set of parameters, a model structure was 
defined on a 200 x 90-cell grid in spherical coordinates, 
{R, Q}. The radial grid was logarithmically spaced from 
the sublimation radius to 1500 AU, with local refinements 
near i?sub, 10 AU (the fixed edge of the inner disk), and 
i?cav to sample the optical depth gradients at those lo- 
cations with higher resolution. The latitude grid was 
linear, with three distinct resolution levels that depend 
on the vertical structure of the model: high near the con- 
densed midplane, coarse at high latitudes near the pole, 
and intermediate in the disk atmosphere. Simulations 
of the stellar irradiation of that model density grid were 
conducted to determine the dust temperature structure, 
using the two-dimensional axisy mmetric JVIonte Carlo ra- 
diative transfer code RADMC (see iDuUemond fc DominikI 
l2004al ). Those calculations explicitly accounted for the 
finite size of the star (i?* ; see Table [3]) to properly treat 
any shadowing effects from the inner rim or cavity wall 
(stellar limb darkening is not considered). For a specified 
viewing geometry (a disk inclination, i, and major axis 
orientation, PA) , a raytracing algorithm was used on the 
results of the radiative transfer calculations to compute 
a synthetic SED and set of 880 /Ltm continuum visibilities 
(sampled at the same spatial frequencies observed with 
the SMA). 

Aided by a series of refined grid searches in parameter- 
space, this initial modeling effort was primarily con- 
ducted manually. A confirmation of fit quality was per- 
formed by comparing the synthetic data for a given pa- 
rameter set with the observed SED and the deprojected, 



elli ptically-averaged 880 /xm visibility profile (for details, 
see [Andrews et al]|2009D . For this model setup and such 
a large sample, parameter estimation with automated 
minimization algorithms would be computationally pro- 
hibitive: the free parameter-space is too large (8 density 
parameters, 2 dust parameters, and 2 viewing geometry 
parameters) and degenerate. IVIoreover, there are non- 
trivial technical obstacles: for example, great care is re- 
quired in defining the model grid for a given parameter 
set, to ensure proper treatment of large optical depth 
gradient s and small-scale fea tures like the rim and wall 
(see also [Andrews et 311120091 ) . Regardless of these chal- 
lenges, good (albeit not unique, best-fit) matches to the 
data were found with a simple and systematic approach. 

In most cases, the viewing geometry was fixed based 
on models of CO spectral images in the literature 
jDutrcv et al. 19981; IPietu et al.l[2007l : llsella et al.l[20T0bl : 
[Lyo ct al. 2011i) or the new J=3-2 data from the SMA 
(for LkHa 330, UX Tau A, and RX J1615-3255; de- 
tails will be presented elsewhere). No such information 
is available for the Ophiuchus disks (SR 21, SR 24 S, 
DoAr 44, and WSB 60): estimates of {i, PA} in those 
cases were inferred from the aspect ratio and orientation 
of the continuum, and are considerably uncertain (see 
lAndrews et al.|[2009ll2010b|) . In practice, the cavity size, 
-Rcav, was determined from the null position in the SIM A 
visibility profile (see §3.5). The dust parameters in the 
cavity, {smax, p}, were set based on the strength of the 
10 /im silicate feature, and the density contrast ((Scav) and 
rim height (/irim) were scaled to match the weak infrared 
excess from the inner disk. The height of the cavity wall 
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(^waii) was adjusted to reproduce the shape and rise of 
the mid-infrared spectrum, and the scale height param- 
eters, {he, ip}, were modified based on the far-infrared 
SED. The characteristic radius, i?c, and density normal- 
ization, were inferred from the shape and amplitude 
of the 880 visibility profile. 

3.5. Parameter Effects and Degeneracies 

To reinforce the intuition used to model these data and 
provide an overview of the individual parameter effects, 
it is helpful to decompose the disk structure into four es- 
sential elements: the inner rim, dust-depleted inner disk, 
cavity wall, and outer disk. Each element has a distinc- 
tive impact on the dust emission spectrum, and together 
they are capable of producing a remarkable diversity of 
observational signatures. To illustrate those signatures, 
we define a fiducial reference model and show separately 
the contributions from each structural element on the 
SED in Figure 2] (fiducial parameters are listed in the 
caption) . Figure [S] explores how the SED and 880 /im 
visibility profile respond to modest changes in the pa- 
rameters of specific interest for transition disks: the rim 
and wall heights (/irim and /iwaii, respectively), the in- 
ner disk mass (set by (5cav), and the cavity size (-Rcav)- 
The effect s of more global structure parameters were dis- 
cussed bv [Andrews et al.l ()2009[ ). 

The emission from the directly-illuminated dust in the 
inner rim dominates the near-infrared SED (see Figure 
4). The luminosity of that hot component is set by the 
size (and optical depth) of the rim, such that a more 
"puffcd-up" rim (larger /irim) produces a brighter excess. 
A sufficiently high (and dense) rim can intercept enough 
starlight to shadow material at larger radii. The atten- 
uated irradiation field that reaches those components 
produces less heating and therefore a fainter infrared 
spectrum (see Figure [5^) . While this shadowing effect 
has been explored for normal disks ([Dullemondl 120001 : 
iDuUemond et al.ll2001[ ). its effects can be pronounced for 
transition disks because much of their infrared spectrum 
is produced b y the directly-illuminated portion of the 
cavity wall (Espail lat et"ani2010( ). lEspaillat etUI ()2011l ) 
have shown that the Spitzer IRS spectra for many transi- 
tion disks exhibit remarkable "see-saw" variability, such 
that an increase in the near-infrared emission is corre- 
lated with a decrease at mid-infrared wavelengths (and 
vice versa). They make a compelling argument that 
this variability is produced by changes in the height of 
the inner rim, as dernonstrated in Figure [5k (see also 
IDuUemond etal|[200l fFlahertv fc Muzeroll3l20lol ) . 

Beyond the rim, the tenuous inner disk emits an in- 
frared spectrum corresponding to the temperature range 
between i?sub and 10 AU. In cases where this material 
is optically thin, the luminosity scales directly with the 
density (parameterized here by <5cav)- However, small 
amounts of dust in the cavity can be optically thick. 
Even for substantial optical depths, the infrared lumi- 
nosity still partly scales with density: for a given h, a 
denser disk has more dust at larger heights above the 
midplane, increasing the height and temperature of the 
infrared photosphere. The result is a brighter (and hot- 
ter) excess spectrum. The same is true for the inner rim, 
and therefore Jcav can also have an impact on shadowing. 
Provided a sufficiently high rim height, and (5cav are 
degenerate (see Figure [5^, b): higher rims can compen- 
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Figure 4. A fiducial model SED decomposed into individual 
structural elements. The ordinates are luminosity densities, = 
iTTiPvF^, in solar units. The model assumes a K3 star (~1Mq, 
2R0) irradiates a disk with a total mass 0.07 Mq, Rc = 75 AU, 
and scale height distribution with -ffioo = 3AU and ?/) = 0.2. The 
densities inside Rcav = 40 AU have been scaled down by a fac- 
tor (5cav ~ 10~^, and the inner rim and wall have been artificially 
puffed up a factor of ^3 over their local H value {H^^^y 0.005 AU, 
or 0.5/?*, and -ffwall ~ 3AU, respectively). The dust in the inner 
disk is the same as in the outer disk atmosphere, with Smax = 1 
and p = 3.5. The solid black curve is the sum of the colored curves 
and the contribution of the stellar photosphere (dashed) . 

sate for lower densities (at least partially). Moreover, the 
inner disk size (fixed here to 10 AU) is also degenerate 
with 5ca.v'- a more extended inner disk (a smaller gap, or 
even no gap) can usually be accomodated by decreasing 
(Scav Other fixed assumptions about the inner disk struc- 
ture can also have comparable effects (e.g., the width of 
the rim, shape of the density profile, local grain proper- 
ties, etc.). The important point to be made is that even 
a small amount of (~^/im-sized) dust inside the cavity can 
have a substantial impact on the SED, depending rather 
sensitively on its detailed spatial distribution and grain 
properties. Therefore, without spatially resolved obser- 
vations, the detailed characteristics of the material inside 
transition disk cavities are strongly model-dependent and 
highly uncertain. 

Despite our ignorance of the detailed structures of the 
inner disk and rim, the cavity size can be measured from 
the SMA data with little ambiguity. Past efforts to in- 
directly infer i?cav relied on estimating the wall temper- 
ature from the shape of the mid-infrared spectrum (e.g., 
D'Alessio et al.i 2005; Calvct ct al. 2005; Espaillat eTaLl 
120071 ). Those models have become increasingly sophis- 
ticated to try and ac count for the irradia tion gradients 
caused by shadowing (jEspaillat et al.|[20Tot ). but still suf- 
fer from the implicit uncertainties in the inner disk prop- 
erties (not to mention the detailed characteristics of the 
wall itself). However, a direct measurement of i?cav is 
possible with a resolved, optically thin tracer like the 
millimeter continuum. The sharp emission contrast be- 
tween the faint dust in the cavity and the bright material 
in the outer disk produces a clear "ring" morphology at 
880 /im. The Fourier transform of that emission profile 
(i.e., the visibility profile) has distinctive nulls at discrete 
spatial frequencies, such that the location of the first null 
(where the real part of the visibilities change sign) is a 
relatively unam biguous measure of the inner edge of the 
ring {Rcav', see 'Hughes et al.l 120071 ). Figure [5}; demon- 
strates the sensitivity of this null position to i?cav; mod- 
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Figure 5. Effects of small modifications to key transition disk parameters on the SED (top) and 880 /^m visibility profile {bottom). The 
black curve is the fiducial model shown in Figure |4l The effects of parameter increases relative to the fiducial model are shown in red, 
decreases in blue: (a) height of the inner rim (/iriin)i a factor of 2; (6) inner disk density (l/(5cav), a factor of 10; (c) cavity size (i?cav), 
±5AU; and (d) cavity wall height {/iwalOi ^ factor of 2. 

ifications on the order of 10% are readily detectable for 
typical noise levels. However, we should note that these 
null positions also have a roughly linear dependance on 
the radial gradient of the surface brightness profile, which 
is set by the product of S and the midplane temperature 
profile. Since the latter is approximately a power-law in 
these models {T^id oc i?""?; typically q « 0.5), the i?cav 
value inferred from a fixed null position is proportional 
to the sum j + q. Since 7 = 1 is fixed in these models, 
we under-estimate the uncertainty on i?cav: allowing 7 
to vary by ±0.5 would contribute an additional ^10-15% 
uncertainty to the cavity sizes. 

In a simple model with an empty cavity, the directly- 
irradiated wall emits a thermal spectrum corresponding 
to the equilibrium temperature at its location (i?cav) 
with a lum inosity proport ional to its vertical extent 
(|Calvet et a l. 2002: D'Alessio et al.ll2005l ). To first order, 
the behavior expected from such a model is qualitatively 
reproduced in Figure [5ji: smaller cavities have hotter 
walls and narrower infrared SED dips. As mentioned 
above, there are slight deviations in the details because 
shadowing can induce vertical temperature gradients in 
the wall. In principle, a high wall can shadow the outer 
disk if the dust is sufficiently settled (low he and/or 
The detailed properties of the cavity wall (e.g., shape, 
thickness) are not well constrained by the SED alone; no 
current observations have sufficient resolution to isolate 
such a (presumably) small feature. 

Although all disk modeling suffers from degenera- 
cies between structure parameters and dust proper- 
ties, it is important to explicitly highlight one key is- 
sue in these transition disk models. The standard 
mass-opacity degener acy f or op tically thin emission 
(e.g., Beckwith & Sar gentl llQQlt [Andrews & Wilfiams 
l2005bl( l and the intricate relationship between irradia- 
tion, h eating, and grain propertie s in the disk atmo- 
spher e p'Alessio et al-IITm [200I I200I iChiang et afl 
120011 ) are important, but will not be reviewed again here. 
The degeneracy of interest is related to how our method 
of inferring densities in the inner disk (i.e., (5cav) depends 



on the dust opacities inside the cavity. Since we lack any 
resolved millimeter signal inside i?cav (although see §4 
regarding LkCa 15), the inner disk densities are deter- 
mined solely from the infrared spectrum. In most cases, 
the relatively strong silicate emission features we observe 
suggest that the cavity is populated by "small" grains 
(i.e., Smax ~ l-lO/xm or less). To be conservative, we 
do not assume that any additional material is present. 
This small dust becomes optically thick even at the rela- 
tively low column densities (^10~^g cm~^, give or take 
an order of magnitude; see Figure |3|) that are present in 
some non-negligible part of the inner disk. The sizes of 
these high optical depth regions are effectively respon- 
sible for setting Sca.v- for a given structure, more mass 
would over-produce the infrared luminosities. In some 
sense, this aspect of the modeling tends to minimize Sca.v 
(or maximize the depletion of dust densities inside i?cav)- 
However, more material could easily be accomodated in 
the inner disk if it has lower infrared opacities (e.g., Smax 
larger than a few /zm; see §4). 

Finally, a few notes about the outer disk parame- 
ters are in order (fo r a more complete discussion, see 
[Andrews et al.l l2009l. It is important to emphasize that 
the nature of transition disk structures makes it difficult 
to constrain radial gradient parameters (7 or For 
smaller disks, the same is true for the characteristic ra- 
dius. Re- The issue is limited dynamic range; only a few 
radial resolution elements can probe the millimeter emis- 
sion in a ring (which is why we chose to fix 7), and the 
far-infrared SED that is less affected by the cavity is only 
sparsely sampled. While we have a reasonable constraint 
on the vertical distribution of dust (he) from the infrared 
luminosities, the current lack of far-infrared data means 
that we do not have a good handle on the scale height 
gradient, ip, for these transition disks with large cavities. 
Linking these models with scattered light observations 
that specialize in constraining the flaring angle (i.e., ^) 
will be beneficial in the future. 

Although they were fixed, the settling parameters {x, 
/} merit a brief discussion. The fraction of large dust 
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grains, /, helps set the height where stellar photons are 
absorbed: smaller / means an increased optical depth 
(more abundant small grains) at a fixed height, and 
therefore a slightly warmer atmosphere that emits more 
infrared radiation. The effects of x, the fractional scale 
height of the large grain population, are subtle. In prac- 
tice, X modifies the vertical temperature gradient at in- 
termediate heights in the disk atmosphere: a larger x 
leads to a cooler temperature in a fixed vertical layer. 
Although potentially important for more direct tracers 
of the disk atmosphere (e.g., scattered light or molecu- 
lar line emission), these effects are indistinguishable from 
adjustments to the scale height parameters, {/ic, ■0}: with 
the data used here. 

3.6. Modeling Synopsis 

The previous sections describe in detail an effort to 
model the SEDs and 880 /^m continuum visibilities for 
this sample of transition disks in a homogeneous, para- 
metric framework. Before presenting those results, we 
first summarize some fundamental aspects of the model- 
ing process: 

• Using two-dimensional, axisymmetric Monte Carlo 
radiative transfer simulations, a parametric dust 
density structure is irradiated by a model star 
to determine an internally consistent temperature 
structure. The results of the simulation are used to 
generate synthetic data products to compare with 
observations. The underlying structure model is 
motivated by a simple viscous accretion disk model, 
with significant dust depletion in a central cavity. 

• Given the relatively faint infrared excesses for tran- 
sition disks, the adopted stellar input spectrum fac- 
tors significantly into inferences about material in- 
side the cavity. 

• Small quantities of (small) dust in the cavity can 
produce large amounts of infrared emission. How- 
ever, the lack of sufficient angular resolution in 
those regions means that the detailed structure of 
that material is unknown; considerable uncertain- 
ties remain for the mass content, spatial structure, 
and dust properties inside the cavity [R < i?cav)- 
Our estimates of density depletion tends to be max- 
imal (Scuv is minimized): more material could re- 
side in the inner disk if it has relatively lower in- 
frared opacities (e.g., grain sizes larger than a few 
/im). 

• Small, local modifications to the vertical distribu- 
tion of material near the sublimation radius (the 
inner rim) and outer edge of the disk cavity (the 
cavity wall) are needed to explain the distinctive 
infrared spectra of the transition disks. Because 
those same features can effectively shadow material 
at larger radii, their properties are at least partially 
entangled with one another as well as more global 
structure parameters. 

• Because the dust-depleted cavities are spatially 
resolved at 880 /xm, their sizes are determined 
directly and with reasonable accuracy (typically 
within ±10% for a fixed 7). 



• Inherent model degeneracies and limited observa- 
tional information make it difficult to quantify sev- 
eral model parameters for an individual disk with 
sufficient confidence. However, because the mod- 
eling analysis was performed in a homogeneous 
framework for the entire sample, intra-sample (and 
qualitative) comparisons of the inferred structures 
arc informative. 

4. RESULTS 

The estimated model parameters are compiled in Table 
01 and the data are compared with the modeling results 
in Figures [6] and [71 Each row in those figures corre- 
sponds to an individual disk, with the panels showing 
(from left to right) images of the data (as in Figure |T]), 
model (synthesized in the same way as the data), and 
residuals (the imaged residual visibilities) , the SED (data 
in black, model in red, stellar input spectrum as dashed 
curve) and the 880 /im visibility profile (data in black, 
model in red). In Tabled we have recast Ec into a total 
disk mass, M^, by integrating the surface density pro- 
file and assuming a gas-to-dust mass ratio of 100:1. The 
scale height normalization (he) is expressed in physical 
units at a fixed radius of 100 AU for easy comparison. 
The inner rim and cavity wall heights are also expressed 
in physical units, the former in terms of the fraction of 
the stellar radius to suggest its potential for shadowing 
material at larger radii. The dust parameters correspond 
only to the material inside the cavity and wall (the dust 
populations in the outer disk were fixed; see §3.2): for 
reference, th eir values in the ISM are Smax = 0.25 liiai 
and p = 3.5 (|D rame fc Leelll984l: iWeingartner fc Drainel 
l2001h . The corresponding model surface density profiles 
are shown together in Figure|8l along with the underlying 
global profiles, Y,g (dashed; see Eq. 1). It is important to 
keep in mind that the surface density gradient was fixed, 
so all profiles have 7 = 1 by default. In practice, this 
means that the inner disk has T, cx l/R, while the outer 
disk is dominated by the exponential taper, E oc 1/e^ 
(when Rc > i?cav, part of the outer disk E profile fiat- 
tens out to vary inversely with radius; e.g., GM Aur). 
The surface densities are also shown in an alternative 
"map" format in Figure [HI which better highlights the 
range of cavity sizes on a linear radius scale. 

The transition disks in this sample have remarkably 
large dust-depleted cavities, with a range of radii (i?cav) 
from 15-73 AU. The sample is likely biased toward large 
cavities, as they tend to produce a more obvious "dip" 
in the infrared SED. However, the extent of that bias is 
unclear, as some disks do not show any obvious signa- 
tures of a cavity in their infrared spectra. The WSB 60 
and MWG 758 disks are good examples; their cavities 
were discovered serendipitously in millimete r-wave im- 
ages (jAndrews et al.ll2009l : ilsella et al.ll2010b[ ). Efforts to 
infer i?cav from the unresolved SED alone can be extraor- 
dinarily challenging (see §3.5). However, there have been 
some notable successes when compared with direct imag- 
ing, including GM A ur ( Calvet e t al. 2005; Hughes et aLl 
|2009ai), SAO 206462 (jBrown et a l. 2007, and this paper) 
and the older disk around TW Hya (Calvet et al.ll2002t 
iHughes et"alll2007f) . While sophisticated SED modeling 
does surprisingly well in many cases — typically agree- 
ing with the i?cav values determined here within ^--^30% 
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Table 4 

Model Parameters 



Name 




Rc 


^^100 




Rcav 


log (5cav 


^rim 




^max 


P 


i 


PA 




[Mq] 


[AU] 


[AU] 




[AU] 




[R.] 


[AU] 


[^tm] 




[°] 


[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


MWC 758 


0.008 


25 


18.3 


0.05 


73 


-5.3 


4.6 


9.5 


0.25 


3.5 


23 


49 


SAO 206462 


0.026 


55 


9.6 


0.15 


46 


-5.2 


1.6 


9.2 


3 


3.5 


12 


64 


LkHa 330 


0.024 


60 


6.5 


0.20 


68 


-5.2 


4.1 


6.8 


1 


2.5 


35 


80 


SR 21 


0.006 


15 


7.6 


0.15 


36 


-5.9 


0.8 


8.2 


10 


3.5 


22 


100 


UX Tau 


0.007 


20 


4.1 


0.35 


25 


-5.7 


0.5 


0.6 


1 


2.5 


35 


176 


SR 24 S 


0.045 


40 


4.5 


0.20 


29 


-5.4 


3.2 


1.5 


0.25 


3.5 


45 


24 


DoAr 44 


0.007 


25 


3.1 


0.02 


30 


-5.1 


3.1 


9.0 


0.5 


3.5 


35 


80 


LkCa 15 


0.055 


85 


2.9 


0.20 


50 


-6.0 


1.0 


7.0 


0.25 


3.5 


49 


241 


RX J1615-3255 


0.128 


115 


3.4 


0.25 


30 


-5.8 


a 


2.0 


0.25 


3.5 


41 


143 


GM Aur 


0.070 


120 


6.2 


0.35 


28 


-5.7 


0.1 


2.8 


0.5 


3.5 


55 


64 


DM Tau 


0.040 


135 


4.2 


0.20 


19 


-4.8 


a 


5.7 


3 


3.5 


35 


155 


WSB 60 


0.028 


30 


10.6 


0.10 


15 


-1.7 


0.4 


0.8 


1 


3.5 


28 


117 



Note. — Col (1): Name of host star. Col. (2): Total mass in the disk (gas+dust, assuming 
a gas-to-dust mass ratio of 100:1). Col. (3): Characteristic radius of the surface density profile. 
Col. (4): Scale height at 100 AU (based on he)- Col. (5): Power-law index of the radial scale height 
distribution. Col. (6): Radius of the dust-depleted disk cavity. Col. (7): Density contrast in the 
cavity, such that iJcav = S/Eg when R < Rcav, where Eg is the surface density profile for a continuous 
disk (see Equation 1 and §3.1). Col. (8): Scale height of the inner rim at the sublimation radius 
(expressed in units of the stellar radius). Col. (9): Scale height of the cavity wall. Col. (10); The 
maximum grain size for the dust population in the cavity and wall. Col. (11): The power-law index 
of the grain size distribution in the cavity and wall, such that n{s) oc . Col. (12): Disk inclination 
(0° is face-on). Col. (11): Position angle of the disk major axis, measured east of north. 
^ We have removed any dust out to a radius of 1 AU for the DM Tau disk and 0.5 AU for the RX 
J1615-3255 disk. Neither case requires a puffed-up inner rim to fit their infrared SED emission (see 
the Appendix). 

(e.g., iBrown et all 120071: lEspaillat et"an 120071 ) - there 
are also some cases where the predicted cavity sizes are 
much less consistent with the SIVIA images. The DlVl Tau 
disk is the mo st extreme example , where the 3AU cav- 
ity inferred by iCalvet et al.l (|2005[ ) is ~6x smaller than 
measured here (i?cav = 19 AU). Such large discrepancies 
are informative; they call for modifications to the as- 
sumed struct ure and /or dust prope rties inside the cavity. 
For example. lEspaillat et al.l (j2011[ ) were able to decrease 
their original estimate of the cavity size in the UX Tau 
A disk by more tha n a factor of 2 (from 71 to 30 AU; 
lEspaillat et "all [2OIOI) using an updated grain composi- 
tion in the cavity wall, both improving their SED-based 
fit quality and providing better agreement with the re- 
solved SIVIA data. These more extreme examples should 
serve as reminders that estimates of cavity sizes from 
SEDs are model-dependent; they should not be consid- 
ered robust until the cavity is spatially resolved. 

The apparent emptiness of the transition disk cavities 
is a striking aspect of the modeling results. Compared 
to the underlying global Eg profiles (dashed curves in 
Figure H]), the surface densities inside these cavities have 
been decreased by ~5-6 orders of magnitude. The WSB 
60 disk is a notable exception, where a more substantial 
amount of dust is required to account for the infrared 
excess around this cool star (the S decrease is a factor 
of ~50, (5cav = 0.02). In most cases, only a few percent 
of a lunar mass of dust (--IO-^^-IO-^IVIq) is required to 
reproduce the observed infrared excesses. As explained 
in §3.5, the dust depletion factors (5cav; and therefore E 
in the inner disk) were effectively determined from the 
infrared spectrum, and therefore are only relevant for 
small (^/im-sized) dust grains in the cavities. Increased 
densities (larger Scav) could easily be accomodated for 



a population of larger dust grains that have lower in- 
frared opacities (i.e., Smax ^ 1/J.m). We can crudely 
estimate the mass of mm-sized grains inside the cavi- 
ties from the 880 fim residual maps in Figures [5] and [71 
Since our models produce no millimeter emission inside 
i?cav, the residual emission near the disk center should 
roughly trace the product of the dust opacity, inner disk 
mass, and a m ean temperature {F^, cx KMca,v{T); e.g., 
lAndrews fc Wi lliams 2005b). In most cases, those resid- 
uals are below the noise floor. Using 3 cr upper limits 
from the residual map centers (see Table [5]) , an 880 //m 
opacity of 3.6 cm^ g~i (see §3.2), and a mean tempera- 
ture of 75 K, the absence of significant residuals suggests 
that there is typically less than ~lMffl (< 10"^]V[q) of 
mm-sized grains in these disk cavitiesO Due to our sensi- 
tivity limits, these 880 /im constraints on the density con- 
trasts are relatively weak {Scav < 0.01-0.1). The LkCa 15 
disk is a remarkable exception, as it shows a significant 
{^6 a- Fi, w lOmJy) and centralized residual peak that 
suggests a small reservoir ('-^IJVIq) of mm-sized grains 
are present inside its cavity. 

Aside from the large dust-depleted cavities, these 
transition disks have surface densi ty profiles similar to 
those i nferred for normal disks (see lAndrews et ani2009l 
l2010bl ). Both populations exhibit the same range of char- 
acteristic scaling radii and disk masses, with Rc — 15- 

^ Note that the constraint from the 880 fim residuals in this 
simple example is technically less stringent than could be made 
from the infrared emission. The 1 /im opacities for the large 
grain population, {smax — 

1mm, p = 3.5}, are only ~50x lower 
than for Smax = 1 (see Figure jSjl, meaning these "mm-sized" 
grains would produce too much infrared radiation before they cre- 
ate 880 /^m residuals (when Afcav J~ 10~*Mq). Nevertheless, the 
example is instructive in pointing out that the 880 fira data are not 
quite sensitive enough to provide very strong limits on Ocav • 
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Figure 6. Comparisons of the disk structure models with observations. The left panels show the SMA continuum image, the model, 
and the imaged residuals, with 3cr contour intervals (see Table [2|. Crosshairs mark the disk center (see §2 and Table [TJ and major axis 
orientation; their relative lengths represent the disk inclination. The right panels show the SEDs and 880 fim visibility profiles, with models 
overlaid in red. The stellar photosphere models are marked with dashed curves. 



135 AU and Md = 0.006-0. 128 Mq inferred here (assum- 
ing a 100:1 gas-to-dust mass ratio). Surprisingly, the 
transition disks are distributed along the same Md-Rc 
correlation that was no ted for normal disks in Ophiuchus 
([Andrews et al.l[2010b( ). with essentially the same scat- 
ter. Of course, disk mass comparisons between these two 
populations are misleading, as the cavities in transition 
disks are substantial mass reservoirs in normal disks. The 
same disk models without dust-depleted cavities (i.e., if 
S = Eg at all radii) would be ^1.3-3x more massive 
than indicated in Table |4] (^20 x for the disk with the 
largest cavity, MWC 758). This implies that the surface 
densities in the outer disk are typically up to a few times 
higher in transition disks compared to their normal disk 
counterparts. 



With a relatively small and potentially heavily-biased 
sample, there is little that can be concluded from any 
trends between the model structure parameters. As 
noted above, the outer parts of the transition disks do 
seem to follow the tentative Md-Rc correl ation identified 
from a sample of normal T Tauri disks ([Andrews et al.l 
2010b). With a brief examination of the figures here 
(with the disks ordered by spectral classification), one 
might notice a possible trend between the stellar type 
and cavity size. Indeed, there is a marginal hint of a pos- 
itive correlation (~2.5(t) between i?cav and T^s (also 
and M*), such that the earlier type stars appear to have 
disks with larger cavities. However, the trend is likely 
an artifact of the endpoints being outliers: the low-mass 
stars host disks with small cavities (DM Tau, WSB 60), 
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Figure 7. Same as Figure|6l The images for DM Tau and WSB 60 are slightly offset to malce room for the inset maps that were synthesized 
with higher angular resolution (see Table [2]l. 



and some of the higher-mass stars have disks with larger 
ones (MWC 758, LkHa 330), but there is a wide range of 
i?cav for intermediate cases. No other trends among the 
disk model parameters, or between them and the stellar 
properties, are found here. A more individualized de- 
scription of the model results for each disk, along with a 
brief comparison to previous models when available, can 
be found in the Appendix. 

5. DISCUSSION 

5.1. Transition Disks and Evolution 

5.1.1. Nomenclature and Our Definition of Transition Disk 

Before integrating this sample of transition disks into 
a larger evolutionary context, it seems necessary to com- 
ment on its placement in the rather disorienting land- 
scape of "evolved" disks in the current literature. Some 



of that confusion is related to nomenclature, and some 
is due to the inherently subjective nature of empirical 
classification. Our definition of a "transition" disk is 
conservative: we mean only those cases with unequivo- 
cal evidence for a sharp drop in the dust optical depths 
inside some radius, either in the form of a dip and 
then dramatic rise in the infrared spectrum over a nar- 
row wavelength interval and/or a clearly resolved cav- 
ity at the disk center. Various other types of disks 
have been associated with the transition phase, includ- 
ing cases without a distinctive dip in the infrared. Such 
disks - sometimes termed "anemic" or "homologously- 
depleted" - have relatively faint infrared excesses and 
often low accretion rates, but do not (yet) exhibit any 
compe Uing observationa l evidence for a dust-depleted 
cavity (jLada et al.l 120061: iNaiita et ai]|2007t iCieza et al.l 
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Figure 8. Model surface density profiles (solid curves), assuming a gas-to-dust mass ratio of 100:1. The underlying global density profiles, 
Tig (see Eq. 1) are overlaid as dashed curves. The E gradient parameter was fixed in the modeling to 7 = 1, meaning that E oc 1/R at 
small radii [R < Rc) and S cx: 1/e^ at large radii (i? > Rc)- Intended as references, the gray squares correspond to the estimated surface 
densities in the Minimum Mass Solar Nebula (MMSN), a representation of the primordial, solar-composition disk densities required to 
account for the current giant planet masses in our solar system (and their uncertainties and feeding zones; Weidenschilling 1977b). Dotted 
lines in the inner disk mark the expected surface densities for 1 earth mass and 1 lunar mass of material (gas + dust), assuming E oc 1/R 
from 0.1-lOAU. Note that the material inside the disk cavities is not resolved with our SMA observations; the projected resolution limits 
correspond to R ^ 20-40 AU (depending on the distance to an individual system). 
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Figure 9. The model surface densities in a map format are shown together with the same color stretch (a E scale- bar is on the right). 
This linear representation of the E profiles provides a clear illustration of the range of cavity sizes, outer disk densities, and sizes (Rc)- 
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[20081 [20I0I: ICurrie'lFaI1[200l ICurrie fc Kemro] 120091) . 
While they may signal an interesting evolution phase 
in their own right, these anemic disks are not bona 
fide transition disks in the sense intended here. In- 
stead, the transition disks (in our adopted meaning) 
are typically identified wit h a set of specialized infrared - 
color criteria d Cicza ot all 120071 120101: [Fang et all 120091: 
iMuzeroUe eral.h2010; .Merrn et al.ll2'010f r However, none 
of these color-selection methods successfully recovers our 
entire sample (with failure rates up to 60%), despite the 
unambiguous confirmations of large disk cavities from 
the SMA (the Men'n et al. "cold disk" criteria perform 
well; a slight expansion of their selection-space would en- 
compass nearly all the disks in this sample). 

These low recovery rates are not surprising, especially 
given the sometimes bright (and complex) emission spec- 
trum generated by a small amount of dust inside the 
disk cavities. Due to the bewildering diversity of possible 
infrared SED behaviors for transition disks, any (unre- 
solved) infrared-based identification scheme is biased and 
will miss good candidates. Unfortunately, our current ig- 
norance of the material and structural properties of the 
inner disk makes it difficult to properly characterize those 
biases. In the end, the only appropriate systematic way 
to find and characterize the transition disk population is 
with a flux- and resolution-limited survey that can di- 
rectly image dust-depleted cavities. Since that type of 
unbiased survey is not yet feasible, we focus on the small 
piece of selection-space that is (at least partially) ad- 
dressed by the specific type of transition disks in this 
sample: namely, those with bright millimeter emission 
and large cavities. 

5.1.2. Frequency of Transition Disks with Large Cavities 

Despite the narrowed scope, transition disks with large 
cavities are surprisingly common when viewed in the 
proper context. Consider the millimeter- wave luminosi- 
ties fr om disks in the nea r by Ta u and Oph associa- 
tions ([Andrews fc Williams! l2005bL [2007at ICieza et al.l 
I2010D . Removing non-detections, sources with likely en- 
velope contamination (including fiat-spectrum objects), 
and disks around early-type (A/B) stars, and scaling 
1.3 mm photometry when necessary by a conservative 
factor 1/ (steeper than the typical slope found by An- 
drews & Williams), we constructed the distribution of 
880 /im luminosities (£mm) shown in Figure [TOl for 91 
disks. Out of that sample, large cavities (i?cav ^ 15 AU) 
have been directly resolved for 9 di sks (10%, includ- 
ing RY Tau; see llsella et al.l l2010al ). That fraction 
is roughly consiste nt with expectations from infrared- 
selection methods ('Muzero lle et al.l 120101) , especially if 
the upper limits for millimeter-faint disks are included. 
But current disk imaging surveys are sensitivity-limited 
and have only probed the brightest disks. In the upper 
half of the I/mm distribution (> 80mJy at 140 pc), the 
transition disks with large cavities comprise a larger frac- 
tion, ~20% (9/46 disks). And in the upper quartile of 
the distribution (> 200 mJy at 140 pc), that fraction is 
slightly higher still (26%; 6/23 disks). These transition 
disk fractions are lower bounds; there could very well be 
more disks with large cavities that have yet to be im- 
aged at sufficient resolution. The maximum fraction of 
this type of transition disk in the upper half (or quartile) 
of the Lmm distribution is ~70%, based on the known 




[Jy] (880 /Ltm at d = 140 pc) 



Figure 10. The distribution of 880 /im luminosities (Lmm), scaled 
to a distance of 140 pc, for disks in the Taurus and Ophiuchus star- 
forming regions (grey), excluding non-detections, sources with po- 
tential envelope contanimation (including flat-spectrum objects), 
and ea rly-type (A/B) stars (data from Andrews & Williams 2005y, 
I2007bl ) . The median and quartiles of the distribution arc shown as 
dashed and dotted lines, respectively. The Lmm distribution of the 
subsample of Tau/Oph millimeter-bright transition disks with di- 
rectly resolved large (i?cav ^ 15 AU) cavities is shown in red. The 
hatched regions correspond to disks that do not have large cavities, 
as determined from millimeter-wave observations with compara- 
ble resolution to this study (Andrews et al. 2009, 2010b). In each 
bin, the fraction of transition disks with large cavities is at least 
(red)/ (grey) and at most Af(hatched)/A(grey). That fraction 
is at least 20% (and at most 70%) for the upper half of the Lmm 
distribution, and at least 26% (at most 70%) for the upper quartile. 
As indicated by the large number of disks that lie in the unhatched 
grey areas of the Lmm distribution, many disks have not yet been 
observed at sufficient resolution to detect a dust-depleted cavity. 
Moreover, we do not yet know the transition disk fraction at low 
Lmm (below the median) or smaller iJcav ■ 

absence of large cavities in s ome of the remaining disks 
([Andrews et al.llM)l l20T0bl ). 

To reiterate those remarkable statistics: at least 1 
in 5, and possibly up to ~two-thirds of the known 
millimeter-bright disks have large dust-depleted cavities 
(with i?cav 15 AU). We do not know how high such 
fractions might be for smaller cavities, or disks with even 
fainter millimeter emission, or young disk populations in 
other environments (if these results for Tau/Oph arc to 
be applied to other clusters, one must keep in mind the 
counting uncertainties that correspond to ±10% errors 
on the transition disk fractions quoted above). Regard- 
less, the implication is that a large fraction of the most 
millimeter-bright - and therefore most massive - disks 
show unequivocal evidence for substantial dust evolu- 
tion in the regions most relevant for planet formation on 
timescales of less than ~l-2Myr (for reference, the me- 
dian and quartile Lmm values correspond to disk masses 
of 0.005 Mq and 0.011 M0, respectively, using the stan- 
dard Andrews & Williams assumptions) . This reinforces 
a similar conjecture made bv iNaiita et al.l (|2007t ) based 
on a group of unresolved Taurus disks, although see 
§5.3.3 regarding their corollary claim of low accretion 
rates. 

5.1.3. Comments on Evolution 

In retrospect, it is striking that these transition disks 
could be such bright millimeter sources in spite of their 
depleted cavities. What were these disks like before they 
evolved, and are such progenitor disks still present in 
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these young associations? If we assume that these disks 
had E = Sg at all radii before depleting their cavities, our 
radiative transfer calculations indicate that Lmm would 
typically be <2x larger (up to 4x for the extreme case, 
MWC 758). Note that the emission increase is not di- 
rectly proportional to the mass increase expected from 
fiUing the cavity with dust (see §4), due to optical depth 
effects. Unfortunately, there are only a few such disks on 
the extreme tail of the Lmm distribution. With such low 
numbers, identifying a parent population for these mas- 
sive transition disks is difficult. It is possible that they 
represent the evolved versions of massive flat-spectrum 
sources (e.g., DG Tau, HL Tau). Or perhaps some of 
them brightened after depleting their cavities: for ex- 
ample, viscous spreading of an initially narrow optically 
thick ring could i ncrease £mm (if the fl ux scales with the 
disk size. Re] see [Andrews et al1l2010 bl. 

The cavity sizes and millimeter luminosities are not 
the only clues available to help characterize transition 
disk evolution, and we would be remiss not to address 
potential hints from the cavity contents. In terms of 
their near-infrared excess luminosities, it is clear that 
not all transition disks are alike: some are faint (in this 
sample, DM Tau, GM Aur, RX J1615-3255, and per- 
haps SR 21) and othe rs are quite bright (e.g., DoAr 44). 
lEspaillat et al.l ()2007D suggested that the latter cases, 
dubbed "pre-transition" disks, may be more representa- 
tive of gaps than cavities. That hypothesis (and accom- 
panying nomenclature) supposes that the transition disks 
with brighter infrared excesses trace an earlier stage of 
the inner disk depletion process. That may very well be 
the case, but the properties of the inner disk are so poorly 
constrained that it is not the only possibility. While 
there is no clear association of infrared excess luminosity 
with the depletion factor, Jcav, we were forced to remove 
dust inside ^1 AU to avoid over-producing the observed 
emission for our DM Tau and RX J1615-3255 models 
(see the Appendix). That is in good agreement with the 
original 'Espail lat et al.l (|2007[ ) conjecture, where the in- 
ner disk densities (or optical depths) and excess levels 
are related. However, the strength of the infrared excess 
is also tied to the vertical distribution of dust near the 
star: larger heights generate more emission (see §3, or 
lEspaillat et al.|[201Cl[ ). Perhaps the transition disks with 
brighter infrared emission also have more "puffed" -up in- 
ner rims. This might suggest that those cases with faint 
excesses have relatively lower inner disk gas densities (or 
gas-to-dust ratios), which do not provide sufficient pres- 
sure to support such extended vertical structures. 

5.2. Potential Dust Depletion Mechanisms 

Up to this point, the discussion has been intentionally 
generic with respect to the physical origins of the dust- 
depleted cavities in transition disks. In a broad sense, 
there are two types of depletion mechanisms: those that 
physically remove or "clear" material from the inner disk 
(decrease S), and those that modify the emission prop- 
erties of the dust (decrease k). Various inner disk evo- 
lutionary processes ha ,ve been assessed in previous stud- 
ies of transition disks (iD'Alessio et al.ll2005HNaiita et al.l 
I2007t lAlexander fc Armitagell2009D . We review several of 
the proposed dust depletion mechanisms in the context 
of the specific type of transition disk in our sample. 



5.2.1. Dispersal by Winds 

Among the proposed clearing mechanisms, mass-loss in 
a photoevaporative wind has excellent potential for clear- 
ing a disk cavity (jClarke et al.l 120011: lAlexander et al.l 
2006 a, b). Energetic photons from the central star can 
ionize the disk surface layers and launch a wind from just 
inside the radius where the thermal energy of t he ionized 
gas exceeds the escape veloc ity (i? « l-10AU: lShu et ahl 
ll99llHollenbach et aLllTOOl . If the inward mass (accre- 
tion) flow rate exceeds the mass-loss rate in this wind, 
changes to the disk structure occur slowly, as expected 
for viscous evolution. But as the accretion flow decays 
over time, it will eventually be surpassed by the mass-loss 
rate in the wind (assu ming the photo ionizing flux from 
the star remains high: lAlexander et a l. 2004, 2005). At 
that point, the accretion flow is redirected into the wind 
and gas is rapidly depleted in a gap near the wind launch 
radius. Isolated from its replenishment source (the outer 
disk), the inner disk material quickly drains onto the star. 
The initially small cavity made in this way exposes the 
disk edge near the wind launch radius to unattenuated 
high-energy irradiation, accelerating the dispersal of ma- 
terial and growth of the cavity. Although this photoevap- 
oration process only directly affects the gas phase of the 
disk (only the very smallest dust grains, s < 0.1 iim, are 
directly swept away in the wind) , lAlexander fc Armitaed 
(2007) demonstrated that a simple prescription for the 
radial drift of solids promotes the simultaneous depletion 
of dust inside the wind-blown cavity. 

In models that consider only UV irradiation, the wind 
plays a minor role in the evo lution of a disk struc- 
ture over most o f its lifetime ([Alexander et al.l l2006bl : 
iGorti et al.l l2009f ). Only late in life, when the disk 
mass and accretion flow have diminished, does photo- 
evaporation dramatically influence the disk properties. 
Observationally, the transition disks produced by UV 
photoevaporation should have empty cavities, low ac- 
cretion rates (M < 10"^" M© yr~^) and low masses 
{Md < 0.005 M(T)), and orbit relativ e ly "old " stars (~few 
Mvr:" lAlexander fc Armitaed I2007L l2009| ). Given the 
large disk masses (Tabled]) and high accretion rates (Ta- 
ble [3]) in this sample, it is clear that UV photoevapo- 
ration is not responsible for creating the large cavities 
in these millimeter-bright tr ansition disks (as noted by 
lAlexander fc Armitagdl2007| ). However, there is growing 
consensus that X-rays drive substantially higher wind 
mass-loss rates, and therefore that photoevaporation 
has a much m ore prominent effect on di s k structures 
at earlv times (lErcolano et al.ll2008l 120091: iDrake et all 
120091: lOwen et al.ll2010D . The earlv onset of X-ray pho- 
toevaporation implies that these disks could still have 
material in their cavities, and may have higher accre- 
tion rate s (up to 10~^M(7) yr~^) and disk masses {Md ~ 
0.01 M©; lOwen et al.ll201lD . To assess the role of X-ray 
photoevaporation, we compare the quiescent X-ray lu- 
minosities {Lx in Table [3]) with the cavity sizes (i?cav), 
disk masses, and accretion rates. In this sample, the 
i?cav and M values for a given Lx are much too large 
to be consistent with their formation via phot oevapora- 
tive mass- loss (as noted by lOwen et al.l 120111) . That is 
not to say that photoevaporation is unable to clear disk 
cavities, as it seems a likely explanation for some evolved 
disks (e.g.. ,Cieza et al..,20G8; Pascucci fc Sterzik.2G09, ) - 
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only that it is insufficient to explain the population of 
massive transition disks with large cavities. 

Some alternative wind sources have been suggested 
as potential contributors to the clearing of disk mate- 
rial. Based on their local MHD simulations of turbu- 
lence driven by the r nagne to-rotational instability (MRI) , 
iSuzuki fc Inutsukal (|2Q09f ) postulated that the disrup- 
tion of large-scale channel flows by magnetic reconnec- 
tion events can drive a substantial wind off the disk sur- 
face. Extrapolated to global evolution models for the 
disk structure, the mass-loss rates in the inner disk from 
such a wind c ould rival those pro duced by X-ray pho- 
toevaporation (jSuzuki et al.ll2010l ). Those initial calcu- 
lations suggest that i?cav ^ 15 AU cavities could be de- 
pleted on Myr timescales, but a more comprehensive set 
of observational predictions that can be used to test the 
feasibility of MRI disk winds as a transition disk clear- 
ing mechanism (and differentiate it from photoevapora- 
tion) is not yet available. Aside from winds launched 
from the disk itself, the wind generated by the central 
star may similarly facilitate disk dispers al by strippin g 
material out of the disk atmosphere (e.g. . iHoredtl 1 1 978l ) . 
Recent models by Matsuyama et al. (200§) indicate that 
the mass-loss rate of disk material from this process is 
quite small for young disks, although perhaps the stel- 
lar wind can act in concert with disk winds to accelerate 
the dispersal of the inner disk. That said, the potential 
synergy of these various wind dissipation models has not 
yet been explored in detail. 

5.2.2. Material Evolution - Particle Growth 

Instead of a lack of mass in the inner disk, the cavities 
we observe could instead be a sign that substantial mate- 
rial evolution has occurred in the inner parts of transition 
disks. Solid particle growth is a natural disk process that 
might be able to explain the transition disk structures 
without requiring drastic density modifications. Because 
large solid particles have low opacities (k), significant 
grain growth can diminish t he c ontinuum emission over 
a wide wavelength in terval (jMiyakc & Nakagawa 1993 
iD'Alessio et al.ir2006[ ). To produce the observed infrared 



SED deficits with a continuous density profile (S = Eg), 
the ~l-10/zm opacities in the inner disk need to be de- 
creased by a factor of roughly ^cav (because the infrared 
emission sets the contrast parameter; see §3.5). The 
880 /im opacity also needs to decrease by a factor of >10- 
100, to account for the ring-like emission observed with 
the SMA. With the dust composition assumed here, such 
opacity spectra are only possible after growth beyond 
Smax = 1 m (preferably with a top-heavy size distribu- 
tion, p < 2.5). Coagulation timescales are expected to 
be shorter in the inner disk, due to the increased densi- 
ties and relative velocities closer to the central star. In 
principle, that would produce a radial variation in par- 
ticle sizes and a decrease in k - and therefore emission - 
toward smaller disk radii. 

Nevertheless, there are challenges to reproducing the 
transition disk cavities with particle growth alone. For 
example, models of the grain size evolutio n in disks pre- 
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visibility nulls) and the shape of the mid-infrared spec- 
trum, particle sizes would need to increase dramatically 
over a narrow radial range (at most a few AU). That 
type of break in k{R) might be reproduced if particle 
growth is significantly enhanced in a large, quiescent 
dead zone with a sharp outer edge. In a dense disk, 
the midplane can be so shield ed from ionizing radiation 
that t he MRI will not opera te ()Gammielll996HSano et afl 
I2OOOI : iFromang et al.ll200^ : without that source of dis- 
ruptive turbulence, particle growth can proceed more ef- 
ficiently inside the dead zone. But there are two prob- 
lems with this scenario. First, the disk atmosphere 
will remain partially ionized and can channel material 
from the outer disk toward th e star in surface accre- 
tion flows (e.g., iGammij 119961 ). Even small amounts 
of fim-sized dust grains carried along in that flow could 
produce enough infrared radiation to wash out the dis- 
tinctive dip in a transition disk SED. And second, the 
cavities are much larger than is expected for dead zones 
(jMatsumura fc Pudritzll2003() . Even in a large dead zone 
without dust in the surface layers, a collisional fragmen- 
tation process might still produce a substantial reser- 
voir of small grains a nd limit growth beyond ^cm -scales 
(|Brauer et all [200I [BiTnstiel et all [2OIOI f20Tll ). New 
models of particle growth under dead zone conditions 
(with radiative transfer) are needed to determine if co- 
agulation and fragmentation in a low-turbulence envi- 
ronment can produce an infrared spectrum and 880 fim 
emission ring consistent with those observed for transi- 
tion disks. 

Molecular emission tracers can provide additional in- 
sight here, as particle growth itself is not expected to sig- 
niflcantly perturb the gas (H2) densities inside the cavity. 
The measured accretion rates (Table [S] aside from SR 
21) unambiguously confirm that neutral gas is fiowing 
through these cavities and onto the stars. Also, some 
of these transition disks emit CO rovibrational spec- 
tra that attest to the p resence of warm molecular gas 
in th e ir inner few AU (iNaiita et all 120031: iSalvk et al"] 
120071 120091 iPontoppidan et all \200&t} . However, new 
studies indicate that those inner regions are deficient 
in other molecule s (e.g ., II9 O: see Naiita et al.l 120101 : 
IPontoppidan et al.|[201Ci[ ). and lDutrev et al.l ()2008l ) claim 
that the CM Aur disk has a i? « 20 AU region of depleted 
CO emission. In some respects, it is clear that transi- 
tion disks have different molecular signatures compared 
to their normal disk counterparts. But given the inher- 
ent challenges in interpreting such measurements, it re- 
mains uncertain whether their cavities also have depleted 
gas densities (other effects like photodissociation, self- 
shielding, and gas-grain chemistry must also be consid- 
ered: e.g..lAikaTya fc Nomura 2006; Jonkhcid et al. 200"^ 
iVasvunin et al.l I2OI 111 . Nevertheless, we can expect that 
these molecular emission lines will be crucial in evaluat- 
ing what role particle growth has in producing the ob- 
served dust structures of transition disks. 

5.2.3. Tidal Interactions with Companions 
A l arge fraction of stars have co mpanions (|Abt fc Levvl 



rather than the abrupt optical depth change near i?cav 
that is observed. To account for the sharp optical depth 
contrast signaled by the 880 /zm emission rings (i.e., the 



ll976HDuquennov fc Mavodll991| ). with stellar multiplic 
ity fractions approachi ng 75% in the diffuse young clus- 
ters n ear the Sun Ce.g.. lWhite fc Ghezll200lHKraus erall 
120111) . Circumstellar material is expected to be strongly 
affected by dynamical interactions with these stellar com- 
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ponents fsee iLin fc Papaloizoul 119931 ) . The interplay of 
resonant and viscous torques can open a gap in a circum- 
system disk, truncate individual disk edges, and poten- 
tially expedite disk dissipation (via accretion or ejection), 
depending on the orbital separation (a), eccentricity (e), 
and mass r atio (q) of the compani o n. Ba sed on the cal- 
culations of lArtvmowicz fc Lubowl ()1994[) , a binary with 
a large separation can retain individual disks similar to 
those around single stars, but with their outer edges trun- 
cated at Ro ~ 0.2-0. 5a. For smaller separations, a cir- 
cumbinary ring structure may also be present, with a 
truncated inner edge at Rcb ^ 1.8-2. 6a. iPichardo et al.l 
(|2008D find similar results, but generalize to arbitrary 
eccentricities such that Ret 1.9(1 -I- e° '^)a. Their cal- 
culations suggest that the inner radius of the circumbi- 
nary ring has only a weak dependence on the mass ra- 
tio, at least for q > 0.1. It is natural to associate the 
880 /im emission rings observed for the transition disks 
with these circumbinary structures. In fact, stellar com- 
panions have been detected inside the cavities inferred 
(indirectly) for a few notable transition disks, includ- 
ing CoKu Tau/4 ([Ireland fc Krausll2008D and the older 
HP 98800 and Hen 3-600 systems lUchida e t al. 2004 
iFurlan et al.ll2007t [Andrews et al.|[2010a|) . 

If the transition disk cavities were created by tidal in- 
teractions with stellar companions, their expected orbital 
separations would be a ~ 0.3-0.6i?cav, for high eccentric- 
ities (e ^ 0.9) down to circular orbits. The range of 
cavity sizes measured here would correspond to compan- 
ion separations of a '--^ 4-40 AU. The multiplicity fraction 
for similar primary star masses in this separation range i s 
-10-20% (|Duauennov fc Mavoi<1991l:lKraus et "aLll2QlTI ) . 
which is similar to, but slightly below, the minimum fre- 
quency of millimeter-bright transition disks with large 
cavities we estimated in §5.1.2. Given how small those 
separations are projected on the sky, ~0. 04-0. 20", di- 
rect searches for compan ions have been limi ted by both 
resolution and contrast. iKraus et al.l ()2011l ) have over- 
come such restrictions in Taurus, using a near-infrared 
aperture-masking interferometry technique coupled with 
the adaptive optics system on the Keck telescope. For 
the four Taurus disks in our sample (as well as RY Tau; 
see llsella et al.|[2010aD . Kraus et al. definitively rule out 
the presence of any companions with q > 0.01 (>0.02 for 
DM Tau) in the sep aration ranges implied by i?cav (see 
also lPott et al.|[2"010D . The same applies to the LkHa 330 
and SR 21 disks (A. L. Kraus 2011, private communica- 
tion). For at least half of the known transition disks with 
large cavities, there are no companions with Ms > 20- 
30Mjiip that could be responsible for the observed disk 
structures. 

Until similar measurements are available for the other 
transition disks, their multiplicity status will remain un- 
certain. However, there are some stark differences be- 
tween the transition disks and a = 4-40 AU multiples. 
Tidal interactions substantially diminish disk masses in 
multi ples with this separatio n rangej Jensen ct al. 1994, 
1996"; Ostcrloh fc Beckwith' [19951: [Andrews fc Williami 
2005b: Ciczactal. 2009). To illustrate that point, 
we compare the millimeter- wave luminosities (imm, see 
§5.1.2) with the projected separatio ns of the known mul- 
tiples in Ophiuch us and Taurus (jBarsonv et al.l 120031 : 
IKraus et al.ll201ll) in Figure [TTb . The transition disks 
are marked with red lines, representing the companion 



separation ranges implied by their cavity sizes. More 
directly. Figure [TTb compares the cumulative Lmm dis- 
tributions of the transition disks and the known 4-40 AU 
multiples (the latter includes upper limits and was con- 
structed using the Kaplan-Meier product limit estimator; 
see Feigelson & Nelson 1985). The two populations have 
remarkably different (or equivalently, Md) distribu- 
tions: the median transition disk is at least 30 x brighter 
(more massive) than the 4-40 AU multiples. Of the 44 
known multiples in that separation range, only one - GG 
Tau - is brighter than 30 mJy at 880 /im. GG Tau is 
exceptionally rare, hosting a circumbinary ring not un- 
like the transition disk structures (it technically meets 
our definition of a transit ion disk in §5.1.1), but with an 
enormous Rcb ~ 180 AU ([Guilloteau et al.lll999f) . In any 
case, the absence of millimeter-bright 4-40 AU multiples 
is striking: it suggests that the specific type of transi- 
tion disk studied here has not suffered from the same 
destructive interaction mechanism. 

Some caution is warranted in over-interpreting the dif- 
ferences between these two populations, as selection bi- 
ases, unknown orbital parameters (e, i), and perhaps 
uncertainties in the t idal truncation models (e.g., see 
iBeust fc DutrevI [20051 ) should also be considered. One 
notable distinction is related to the mass ratio. There 
are stringent limits on q (> 0.01-0.02) for at least half of 
the transition disk sample, but the 4-40 AU multiples are 
prefe rentially equal-mass (median q w 0.7; IKraus et al.l 
[20T1I ). Such high mass ratios could impart more severe 
tidal perturbations to the circumstellar environment at 
an early evolutionary stage. Since that seems unlikely 
for the transition disks, it is worthwhile to explore the 
effects of companions with much lower q values. Such 
low-mass companions that form in the disk around the 
primary - including giant planets and brown dwarfs - can 
modify the disk structure without dispersing too much 
of the total mass (oc Lmm). 

Resonant torques induced by a low-mass compan- 
ion set up spiral wave s that repel l ocal disk material 
(Goldrcic h fc Tremaind 119781 Il980t [LmjLPapaloizo^ 
[l979al) bl). Above a mass ratio threshold, q > 
ggap ft! 40Q(j?/a)^ (where a is the viscosity coeffi- 
cient; iShakura fc Sunvaey 1 19731) . these tidal interactions 
open a gap m the disk (ILin fc Papaloizoul 119861 119931: 
iBrvden et allll999l: iCrida et al.ll2006[) . In this scenario, 
the companion is located closer to the outer edge of the 
cavity than for the high-q case. A typical separation 
for a circular orbit might be a — 0.8i?cav, although the 
exact location depends on some complicated issues that 
influe nce the gap width ([Winters et al.|[2003t ICrida et al.l 
[2OO6I) . For our model structures and a = 0.01, the tran- 
sition disks in this sample would need companions with 
a — 10-60 AU and mass ratios at least ggap ~ 0.001- 
0.03 to open gaps near their cavity edges. Using the A/, 
values in Table [3J those limits translate to companion 
masses Ms > 1-70 Mjup (the high end corresponding to 
larger M,; note that (jgap and Ms scale linearly with the 
assumed a). For the relevant cases, the mass ra tios im- 
plied from current infra red contrast limits (Pott et al.l 
I2OIOI : IKraus et al.l 120111 ) exceed those needed to open 
gaps. Although reassuring, there is little to be learned 
there, given the uncertainties on the disk viscosities (a) 
and the luminosities of such low- mass companions (e.g., 
,Marlev et al...2007. ). 
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Figure 11. {left) The 880 fim lumin osity (Lmm; [Andrews He Williamsl [2005bl. |200 7b|) as a function of the projected separation (a; 
iBa rsonv ct al. 2003; Kraus et al.l 120111 ) for multiple star systems in the Taurus and Ophiuchus star- forming regions. Gray points with 
arrows signify 3o" upper limits on Lmm. The red lines mark the stellar companion separation s that would be expected to clear the transi- 
tion disk cavities: their vertical extents correspond to the range of eccentricities e ~ 0.0-0.9 IIArtvmowicz &: Lubowlll994l : IPichardo et al.l 
[2009) . (right) A direct comparison of the cumulative Lmm distributions for the transition disks (red) and the multiples with projected 
separations of 4-40 AU (black; the latter was constructed with the Kaplan- Meier estimator to include 3 a upper limits) . The two populations 
have extremely different millimeter-wave luminosity (disk mass) distributions: transition disks are typically a factor of '~30x brighter. 



An individual gap is too narrow to be responsible for 
the observed properties of the transition disks. How- 
ever, the disk-companion interactions can also regulate 
the densities in the inner disk, interior to the compan- 
ion orbit. As in the case for a photoevaporative flow 
(§5.3.1), the gap effectively isolates the inner disk from 
its replenishment reservoir at larger radii. Without a 
viscous flow of material inwards, the inner disk ma- 
terial would drain onto the star on a relatively short 
timescale, producing a large and empty cavity. Unlike 
the photoevaporation mechanism, some material from 
the outer disk is expected to seep through a gap cre- 
ated by a low-mass companion in accretion streams, 
but only at decreased flow rates compared to a nor- 
mal disk (e.g ., Artymowicz fc Lubow 1996; Kiev 1999; 
iLubow et allTl999h . E^bow fc D'Angeld (2006) estimate 
that the mass flow across such a gap is diminished by 
~75-90%, depending on q: a more massive compan- 
ion reduces the flow across the gap, leading to a lower- 
density inner disk. A quantitative association between 
the companion mass, the accretion flow, and the deple- 
tion of the inner disk is not yet clear in simulations. 
Some have suggested a companion with q sa 0.001 is 
sufflcient to deplete a disk cavity like those observed 
for the transition disks (iRice et al.ll2003at iQuillen et al.l 
120041 ; [\%rniere et al.ll2006( ). while others have argued that 
h igher masses are req uired (iCrida & MorbidcUi 2007). 

INaiita et al.l ()2007l ) suggested that transition disks ex- 
hibit low accretion rates compared to their normal coun- 
terparts with similar disk masses, consistent with a large 
reduction in the mass flow acr oss a gap produced by 
a low-mass companion (see also [Alexander fc Armitagd 
120071 12009[ ). To verify that result, we constructed the 
distribution of M for normal disks in Taurus and Ophi- 
uchus that lie in the bright half of the Lmm distribution 
and compared it with the accretion rates for the transi- 
tion disks with large cavities in those regions (including 
RY Tau; see Table [3]). Figure [12] shows the cumulative 
distributions of M for both of these populations. A two- 
sided Kolmogorov-Smirnov test (or its survival analysis 



eq uivalents if we inc lude the M upper limit for SR 21; 
see llsobe et aLlfTgSGt ) indicates that the probability these 
M populations are drawn from the same parent distri- 
bution is low (<l-2%), with the transition disks having 
about ~3x lower accretion rates on average. That rela- 
tive M decrease (^^70%) is similar to the expected reduc- 
tion in flow rates across a gap opened by a q » 0.0 01 (i.e., 
Ms w IMjup) companion (ILubow fc D'Angelol 12006). 
However, the samples are small, incomplete, and likely 
biased. Moreover, the individual accretion rates have 
large uncertainties and are estimated from different diag- 
nostics. So, although this preliminary comparison (and 
that bv INaiita et al.ll2007[ ) provides a tantalizing hint of 
reduced mass flow rates for transition disks, we caution 
that a more rigorous and homogeneous set of compar- 
isons is warranted before any conclusions are drawn. 

Regardless of those comparisons, the M values for the 
transition disks are still high enough to be seemingly in- 
consistent with the low dust densities we infer inside their 
cavities. One way to reconcile those measurements is to 
include additional low-mass compa nions inside the ob- 
served disk cavity. ,Zhu et al.l (j2011t ) have suggested that 
multi-planet systems can dynamically clear large cavi- 
ties while also rapidly transporting gas across the large 
distances between the remnant outer disk and the star. 
Alternatively, a single low-mass companion might be re- 
sponsible for creating a cavity that is depleted in dust, 
but not in gas. Any companion that can open a gap will 
generate wakes just outside its orbit, creating a pressure 
bum p that may act like a dust filter on the flow across the 
gap |Rice et al.ll2006l; iPaardekoo per fc Melie^l2006l ). If 
that filter removes a substantial fraction of the dust from 
the flow, the inner disk can be preferentially depleted of 
dust (producing the thermal deficits we observe) while 
still maintaining high gas accretion rates. 

The tidal interactions between disk material and low- 
mass companions offer a compelling explanation for the 
transition disk structures. However, the observations are 
only starting to provide some tangible constraints that 
can be incorporated into the numerical simulations. As 
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at wide separ ations from the pr i mary (lBossll2006t iBolev 
2009; RafikOT| [200llClarkell2009l:lDodson-Robinson et al 
200S ). Of course, it is not clear that the transition disks 
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Figure 12. The cumulative distribution of the accretion rates (M) 
for the normal disks in Tau and Oph that reside in the bright half 
of the Lmm distribution (black, or blue excluding the spectroscopic 
binary DQ Tau), compared with that for the transition disks with 
large cavities in those same regions (red; including RY Tau). The 
latter was constructed with the Kaplan-Meier estimator to include 
the upper limit on M for the SR 21 disk. We do not include the 
candidate transition disk RX J1633-2422 in either distribu tion, al- 
though it has a very low M estimate (see Cieza ct al. 2010^ . Accre- 
tion r ates w ere compiled from Hartmann ct al. ( 1998), Nat ta et al.l 
jiaoOej ). and lNaiitaeFal.! 42007). Of the 46 millimeter- bright disks 
considered, 7 (primarily in Oph) have no M measurements in the 
literature. In this preliminary comparison, the transition disks 
typically have accretion rates ~3x (0.5 dex) lower than the normal 
disks. However, it is difficult to reach any definitive conclusions 
from such a comparison: the two samples are small and likely in- 
complete (and biased), and the M values have large uncertainties 
and were determined in different ways. 

the data improve - both for the disk structures and com- 
panion searches - we should be able to make a more 
nuanced assessment of whether the mechanism respon- 
sible for transition disk cavities is tied directly to the 
formation and early evolution of low-mass brown dwarfs 
and/or exoplanets. 

5.3. Implications for Planet Formation 

To summarize the discussion in §5.2, we have argued 
that the dust-depleted cavities in the specific type of 
transition disks studied here - those with large Lmm 
and i?cav - were not likely produced by photoevapora- 
tive winds or tidal truncation by stellar-mass compan- 
ions. The possibility that the transition disk cavities are 
the manifestation of substantial particle growth in large, 
low-ionization dead zones remains, but requires a more 
detailed modeling effort to determine if some significant 
challenges could be mitigated. While these mechanisms 
all play non-trivial roles in disk evolution, the weight 
of the observational evidence for the kind of transition 
disks studied here suggests that the cavities are likely 
the byproduct of dynamical interactions with long-period 
low-mass companions - either brown dwarfs or giant plan- 
ets. If that is the case, it is of great interest to understand 
how those companions can form at (or be moved to) such 
large separations from their host stars (a ^ 10-60 AU) on 
such a relatively short timescale (at most ^1-2 Myr). 

Giant planets or low-mass brown dwarf companions 
could be formed relatively rapidly from a fragm enta- 
tion p rocess in a gravit ationally unstable disk (e.g.. lBossl 
119971 119981: iRice et al.ll2003bD . Moreover, this disk in- 
stability mechanism seems to naturally favor formation 



initially had sufficiently high densities (and efficient cool- 
ing) to be unstable, and in most cases (8/12 in this sam- 
ple) the orbits of any companions would be small enough 
(a < 35 AU) that the prospe cts for fragmentation are 
debatable (but see lBossll2009[ ). If the companions are 
not too massive, the more standard core accretion mech- 
anism for giant planet formation may be applicable. In 
that scenario, the coUisional agglomeration of solids must 
produce a massive core to accrete a massive atmosphere 
before the gas disk disperses (e.g., iPoUack et al.l 119961 : 
iHubickvi et al.l I2005D . Traditionally, that process has 
been considered too inefficient to make long-period plan- 
ets (g > 10 AU ) in situ within the disk lifetime. However, 
IRafikov' (20111) has argued that an improved treatment 
of planetesimal accretion shows expedited core growth 
at low surface densities, thereby enabling the capture of 
a massive atmosphere in the outer disk. His estimate 
of the >0.1g cm~^ of solids needed for that process is 
strikingly similar to (or below) the E values we derive 
near i?cav (see Figure [5]). Others have suggested that 
long-period companions could be produced relatively eas- 
ily by the core accretion mechnism in the inner disk, 
and then s ubsequently moved outward via planet-p lanet 
scattering fSch arf fc Menoull2009HVer as et al."2009') or a 
resonant migration process (|Masset fc SncUgrove 20011 
ICrida et al.il20'09l) . 

5.3.1. Future Prospects 

At this stage, it is not clear that the properties of tran- 
sition disks can help differentiate between these planet 
formation mechanisms - even if it can be confirmed that 
their structures are produced by planetary-mass compan- 
ions. However, a future census of companions embedded 
in such disks that determines the joint distribution of 
{a, Ms} would help guide that debate. In that context, 
the highest priority should be placed on continuing the 
high-contrast and high-resolution infrared searches for 
low-mass companions inside the known transition disk 
cavities, at the very least to rule out the presence of 
t heir more massive ( stellar) counterparts (along the lines 
of lKraus et al.|[20lll) . The identification of any young ex- 
oplanet candidates still embedded in their natal disk ma- 
terial would be a landmark discovery for both planet for- 
mation and disk evolution. Moreover, those companions 
would provide youthful touchstones that could help fa- 
cilitate a better understanding of their much older coun- 
terparts, including the Fomalhaut and HR 8799 planets 
that orbit j ust inside their own remnant debris disk rings 
(IKalas et a l. 2008; Marois ct al. 2008, 2010). 

While those companion searches continue, there are 
two issues related to the planet formation process that 
can be better addressed with observations of the transi- 
tion disks themselves. First, the Atacama Large Mil- 
limeter Array (ALMA) project starts operations soon 
and will routinely achieve extraordinary millimeter con- 
tinuum sensitivity with sub-arcsecond resolution. Those 
capabilities will permit the first large, unbiased, and di- 
rect searches for large (i?cav > 10 AU) dust-depleted cav- 
ities in nearby disks, and should provide a more defini- 
tive assessment of transition disk properties across the 
entire disk mass distribution (see §5.1.2). These imag- 
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ing surveys will provide a first look at the fraction of 
disks that have formed long-period giant planets within 
~l-2Myr, and can be compared with direct imaging ex- 
oplanet searches around older stars in the solar neigh- 
borhood. As those surveys are building up statistics, a 
significant eflFort should be invested in developing our un- 
derstanding of the structure and contents of the material 
inside the transition disk cavities. With a focus on indi- 
vidual sources. ALMA will be sensitive enough to provide 
stringent constraints on the depletion of mm-sized grains 
and molecular gas in these cavities, and could potentially 
detect faint emission signals that can be related to the in- 
ner disk mass. The residual 880 fiia signal in the LkCa 15 
disk cavity is a tantalizing hint of the prospects for such 
studies. Coupled with estimates of the cavity size and 
accretion rate, that information can be used in concert 
with numerical simulations to help infer the mass and or- 
bital period of the planetary companion perturbing the 
disk structure. While the present study is just one of the 
first steps, ultimately high angular resolution radio map- 
ping of disk structures can serve as a new tool to guide 
searches for young exoplanets, as well as to characterize 
some of their fundamental properties. 

6. SUMMARY 

We have presented a sample of new and archival SMA 
observations of the 880 fj,m continuum emission from pro- 
toplanetary transition disks at resolution scales down to 
--O'.'S (20-40 AU). Using two-dimensional Monte Carlo 
radiative transfer calculations and a parametric formu- 
lation for the disk densities, we simultaneously modeled 
these SMA visibilities and the broadband SEDs for each 
disk in a homogeneous framework. The data and those 
modeling results were then synthesized to discuss the 
generic properties of transition disks from a primarily 
millimeter-wave perspective, as well as to assess the roles 
of various potential disk evolution mechanisms. Our key 
results include: 

1. In all cases, the SMA data resolve large (i?cav = 15- 
73 AU), dust-depleted cavities at the disk centers. 
Unlike estimates based on models of the unre- 
solved infrared spectra, these direct imaging mea- 
surements constrain the cavity sizes with little am- 
biguity: i?cav uncertainties are typically only ~10% 
(or up to 15-20% if a range of surface density pro- 
files are considered). 

2. Some dust is still present inside these disk cavi- 
ties, although its mass and spatial distribution are 
poorly constrained. With respect to continuous 
density models, we estimate from the SMA data 
that the densities of mm-sized particles are >10- 
100 X lower inside the cavities. From the infrared 
SED, we find that ^m-sized particles are depleted 
by larger factors, typically 4-6 orders of magnitude. 
Despite these low densities, the dust inside the cav- 
ities can have profound effects on the infrared spec- 
trum and the structures we infer at larger radii. 

3. Outside the cavities, the transition disk structures 
appear to be similar to disks with continuous den- 
sity distributions: they have the same ranges of 
characteristic radii {Rc) and total masses (M^). 



In their mass reservoirs at large radii, the surface 
densities for the transition disks are slightly higher 
than normal disks with similar millimeter luminosi- 
ties. 

4. Our sample is biased toward millimeter-bright 
disks with large (i?cav ^ 15 AU) dust-depleted cav- 
ities. However, we find that this specific type of 
transition disk is surprisingly common among the 
Taurus/ Ophiuchus disks with large millimeter lu- 
minosities (or disk masses): they comprise at least 
20% of the disks in the upper half of the imm (or 
Mfi) distribution, and >25% in the upper quar- 
tile. Imaging siirveys of other millimeter-bright 
disks rule out large-cavity transition disk fractions 
in excess of 70%, but there is no such information 
available for smaller cavities or fainter disks. The 
transition disk fractions derived from direct images 
of cavities are higher than those estimated from in- 
frared colors because the latter technique system- 
atically misses cases with excess emission from a 
small amount of dust inside the cavity. 

5. The large cavity sizes, relatively high accretion 
rates, and low X-ray luminosities suggest that the 
transition disk cavities in this sample were not 
cleared by an internal photoevaporation process (as 
was already argued by the developers of those dis- 
persal models). 

6. If the diminished signatures of dust grains in the 
cavities are to be explained by particle growth 
alone, the transition disk structures need to con- 
tain very large MRI-inactive dead zones, have ex- 
perienced growth to at least meter-sized solids, and 
harbor only very limited reservoirs of small dust 
grains in their inner disks (from collision fragments 
or dust coupled to a layered accretion flow). While 
those stringent conditions make this mechanism 
for cavity formation seem challenging, a real as- 
sessment of the role of particle growth in produc- 
ing transition disk structures cannot be made until 
more detailed models are explored. 

7. In principle, dynamical interactions with stellar 
companions on a « 4-40 AU orbits would be capa- 
ble of opening the disk cavities we measure. How- 
ever, recent infrared searches definitively rule out 
any such companions with mass ratios q > 0.01 
for half of this sample. Moreover, the transition 
disks studied here have substantially higher (typi- 
cally >30x) millimeter- wave luminosities (or disk 
masses) than the disks in multiple star systems 
with 4-40 AU separations, suggesting that they did 
not suffer the same level of mass loss from tidal 
perturbations. 

8. We find the best overall agreement with the data 
if these transition disk cavities are opened by tidal 
interactions with low mass ratio (q > 0.001) com- 
panions - either low-mass brown dwarfs or giant 
planets. Such a companion can open a gap and 
substantially deplete material interior to its orbit 
without disrupting the outer disk. We demonstrate 
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that the typical accretion rates for these transi- 
tion disks are decreased (by ^70%) compared to 
their high-Lmm normal disk counterparts, commen- 
surate with the reduced mass flow that is expected 
to cross the orbit of a low-mass companion (how- 
ever, a more in-depth analysis with larger samples 
is warranted). In the near future, radio imaging 
studies of transition disk structures could provide 
a new and unique method for identifying and char- 
acterizing very young (^^IMyr) low-mass compan- 
ions - including giant planets on long-period orbits 
- still embedded in their natal disk material. 
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APPENDIX 



COMMENTS ON INDIVIDUAL SOURCES 

This appendix highlights some aspects of the modeling results for individual disks, and compares with previous 
efforts to characterize their disk cavities whenever available. Table [5] is a compilation of the original sources used to 
construct the SEDs displayed in Figures [6] and [71 

MWC 758 This Herbig Ae star is relatively isolated, located east of the main Taurus- Auriga clouds. Without 

a firm cluster associat ion, its distance is uncertain but assumed here to be the nominal Hipparcos va, lue, d = 200 pc 
(jPerrvman et al.lll997[ ). We adopt the A8 spectral classification advocated bv lBeskrovnava et al.l ()1999[ ). but note that 
others have utilized an earlie r type (A5;lvan den A nckcr et al. 1998). The SMA observations from the E configuration 
are the same as presented bv llsella et al.l (|2010b[ ). although they have been re-calibrated with an amplitude scale that 
amounted to a ~30% decrease in flux levels. That modification was required to reconcile those data and our new 
SMA observations with short antenna spacings (C configuration); the combined data are now in excellent agreement 
with single-dish photometry measurements. The 880 /im visibility null calls for the largest cavity in the sample, 
^cav = 73 AU, despite the bright infrared excess. The standard infrared signatures of a transition disk are heavily 
muted, washed out by only a small amount of dust inside the cavity. Although that dust has densities that are 
comparable to the other disks in the sample (and only modest rim and cavity heights), it produces a significantly 
brighter infrared spectrum for two reasons: (1) the central star is extraordinarily hot and luminous, substantially 
heating more dust, and (2) the scale heights are large, increasing the irradiated surface area and therefore the emitted 
fl ux from the disk atm osphere. 

Ilsella et al.l ()2010bf ) also modeled the MWC 758 disk structure with some of these SMA data (and supplementary 
3mm observations), choosing not to modify the accretion disk model to include a cavity (i.e., S = Sg at all radii in 
their work) but to allow the gradient of the viscosity profile (7) to vary freely. They can reproduce the data with a large 
negative viscosity gradient, 7 « —4, which results in a surface density profile that increases like inside ~70AU, 
peaks near 100 AU, and then drops off like 1/ exp(it!^) at larger radii. Roughly speaking, that density distribution 
is similar to the narrow ring employed here (see Figure [5]) , even peaking at essentially the same location and surface 
density. Not unexpectedly, the Isella et al. model requires an additional dust component near the star to account 
for the infrared excess (which their i?** density profile alone vastly under-produces), similar to the inner disk material 
utilized here. Finally, the sharp density cutoff at large radii that is inherent to models with negative 7 values means 
that the Isella et al. model has effectively no mass outside R ~ 150 AU. On the contrary, our model has gas densities 
in excess of 10~^ g cm~^ out to i? ss 400 AU, despite the small characteristic radius (because 7 = 1). We expect that 
this difference in the effective disk size would help explain the strong positive residuals at large radii noted by Isella 
et al. when applying their dust model to the CO emission. 



SAO 206462 Commonly referred to as HD 135344 B, this star is also isolated, perh aps associated with the outskirts 

of the Sco-OB2 asso ciation. We adopt th e distance estimate of lvan Boekel et all (I2005D. d = 140 pc, a nd the F4 spectral 
type determined bv lDunkin et al.l (|1997D . The SMA data were first presented bv lBrown et al.l ()2009D . and only include 
observations in the V configuration with few short baselines. Our models imply a large cavity, i?cav = 46 AU, with a 
small amount of dust near the star. The super-heated inner rim near the dust sublimation radius is responsible for the 
substantial near-infrared bump in the SED. To minimize the strength of the 10 /im silicate feature, comparatively large 
(smax = 3 fim) dust grains were used to populate the inner disk. Our model tends to over-predict the flux densities 
at wavelengths longward of 1 mm, a not uncommon feature due to our decision to flx the grain size distribution in 
the midplane of the outer disk (see §3.2). The issue might also be alleviated if new data with short antenna spacings 
become available, enabling a more robust disk mass estimate. The SMA image of the SAO 206462 disk is notably 
asymmetric, with enhanced emission to the southeast and a pronounced deficit to the north. But the imaged residuals 
indicate that those asymmetries are only marginally significant: given the challenges of observing such a southern 
target, those feature s are not inconsist ent with noise (but note that a similar pattern is present in the scattered light 
ir aages present e d by Grady et al.l 120091 ) . 

iBrown et al.l (12007) al so used RADMC calculations to reproduce the SED alone, and later combined with these SMA 
data (jBrown et all 120091 ). but with a different density model (a power-law E profile, truncated at the edges). They 
estimated roughly the same cavity radius, 39-45 AU, consistent within the expected ^10% mutual uncertainties. 
Although their lower i?cav value was determined from the SMA data, the difference with the cavity size derived here 
is likely due to small offse ts in the adopted disk centroids and viewing geometries. Regarding the latter, we utilize the 
inclination de termined by Lvo et al.l (|2011[ ) from CO spectral line images, in good agreement with the scattered light 
constraints of lGradv et al.l (|2009[ ). Those scattered light images suggest a substantially settled vertical distribution of 
dust with an effective ip = —0.3, unlike the more common ip — 0.15 used here. Part of that difference is related to 
shadowing from the inner rim/disk and cavity wall, which can effectively steepen the scattered light surface brightness 
profile. But there is so little radial dynamic range available from this disk that it should be possible to find anti-flaring 
models that are also commensurate with the SMA data. 



LkHa 330 Associated with t he Perseus mole c ular c louds {d w 250 pc; see lEnoch et al.|[2006f) . this young star was 

classified as spectral type G3 bv iCohen fc Kuliil (|1979[ ). The SMA data in the V configuration were first presented 
bv I Brown et al., (,2008 ): new SMA data in the C configuration are also included here. Note that there is considerable 
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disagreement between the IRAS and Spitzer data in the far-infrared, perhaps due to molecular cloud or cirr us con- 
tamination in the former. When re-calibrating the long-baseline SMA data used by I Brown eFall ([2008l[2009l) . it was 
noted that the observations of 2006 November 11 had such poor phase stability that they could not be used. With 
their removal and the addition of new data with short antenna spacings, a much more symmetric and sensitive map 
of the 880 fim continuum emission is available. 

We measure a large cavity, i?cav = 68 AU, ^-^35-60% larger than estimated by iBrown et al.l ()2008l l2009f ) from (some 
of) the same data. Most of that discrepancy is caused by an inaccurate estimate of the 880 /zm visibility null location, 
due to the poor quality of one of the SMA tracks in the V configuration. Using the CO J=3— 2 spectral line images 
from the new data in the C configuration, we also use an updated viewing geometry (with a significant difference in 
the major axis PA; the full CO datasets will be presented elsewhere). The faint 10 /im silicate feature was difficult to 
reproduce with the typical grain size index {p = 3.5), but we found some success using a more top-heavy distribution 
(p = 2.5) for the dust inside the cavity. To reproduce the relatively bright infrared excess with that grain population, 
a substantial increase in the inner rim height was required (the dust size distribution parameters and rim height are 
strongly degenerate; see §3.5). 

SR 21 This young star with spectral t ype G3 is still embedded in the L1688 da rk cloud in the nearby Ophi- 

uchus star-forming region, d ~ 125 pc (e.g.. iLombardi et al.ll20d8L iLoinard et al.ll2008| ). The extinction estimate of 
iPrato et"al1 ()2003l ). Ay ~ 9, is too large to be appr opriate for both t he optical photometry of iVrba et al.l ([l993) and 
the near-infrared magnitudes in the 2MASS catalog (jCutri et al .''2003) . Assuming variability is not responsi ble for that 
discrepancy, we adopt a lower extinction {Ay = 6.3) and the refore a substan t ially l ower L„ compared to .Brown et al.l 
(|2007l I2009f ). The SMA data are the same as described bv [Andrews et al.l ()2009D . although the 880 /im continuum 
map has been synthesized with slightly improved resolution. The 880 fim continuum emission is distributed in a nearly 
face-on ring, with i?cav = 36 AU. The narrow ring width (Re is only 15 AU) makes it difficult to determine the viewing 
geometry, and there is no uncontaminated CO emission or scattered light image available to facilitate an independent 
measurement. The lack of a clear silicate feature in the mid-infrared leads us to adopt a large Smax (10 /im) for the 
dust inside the cavity. Despite th e different modehng se tups, the disk structure parameters inf erred here are essentially 
the same as those determined bv [Andrews et al.l ()2OO90 . With (some of) the same SMA data. IBrown et al.l (20091 also 
find a comparable cavity size (33 AU) . We find positive 3 a residuals for this model, suggesting excess emission to the 
southeast and near the stellar position. Those discrepancies lead to a poor prediction of the secondary visibility null 
(see Figure [6|) , which may indicate that there is a small amount of ^mm-sized dust emitting from inside the cavity. 

UX Tau A The primary star of a wide multiple system in Taurus (d » 140 pc; se e iTorres et al.ll2009l). UX Tau A 

is the only component with evidence for remnant disk material (e.g.. iWhite fc Ghezll2Q01l : iMcCabe 'eralll2006h . No 
880 /jm emission above 4.5 mJy (3ct) is detected around either the UX Tau B or C components (in agreement with 
iJensen fc Akesonll2003l ) . The SMA data presented here are the first resolved observations of the UX Tau A disk. This 
emission is quite compact: we use a small characteristic radius {Rc = 20 AU), presumably produced by the truncation 
of the outer disk by UX Tau C (located ^2'.'5 to the west). Unlike the compact continuum, we find a bright, resolved 
CO J=3— 2 disk orbiting UX Tau A, providing the first direct estimates of the disk inclination and position angle. As 
discussed in §5.1, the 880 /tm image shows a res olved cavity that is substan tially smaller, -Rcav = 25 AU, than initial 
estimates from the SED alone (56 and 71 AU bv lEspaillat et al.ll2007l 120101 respectively). Like LkHQ;330, we assume 
a top-heavy dust size distribution {p = 2.5) in the cavity to avoid over-producing the 10 /im silicate feature. 

SR 2 4 S Embedded in the L1688 dark cloud in Ophiuchus, SR 24 is a hierarchical triple system with a close binary 

(0'.'2: lSimon et al.lll995l) located ^5" north of a spectral type K2 primary (SR 24 S). The SR 24 N binary exhibits some 
thermal excess presumably associated with a low-mass circumbinary disk, although it is not detected at millimeter 
wavelengths (despite evidence for associated CO line emission; Andrews & Williams 2005a). Due to potential confusion 
with the SR 24 N component, we did not use a Spitzer IRS observation of this system (although one apparently exists). 
The 880 /im image of the SR 24 S disk shows an elongated ring with a significant brightness asymmetry along the major 
axis (brighter to the north). We me asure a moderate cavity size, i?cav = 29 AU, essentially the same as our previous 
effort (32 AU: [Andrews et al.ll2010bf ). Naturally, the model does not reproduce the observed asymmetry, instead over- 
predicting the emission to the southwest (although marginally, based on the imaged residuals). In general, the fit 
quality here is superior to our initial work, thanks to the modification of the inner rim/disk structure. However, 
without a detailed infrared spectrum the dust properties and structure of the cavity are uncertain. Unfortunately, the 
asymmetry of the narrow emission ring makes a clear viewing geometry estimate difficult, and the CO emission from 
SR 24 S is contaminated by local cloud material. Substantial improvements in the structure constraints are feasible 
with a more complete SED, if the disk inclination can be measured independently (e.g., with a tracer like ^'^CO). 

DoAr 44 Sometimes referred to as ROXs 44, this young star is located in the densely -populated L16 8 8 dar k 

cloud in Ophiuchus. A preliminary analysis of the DoAr 44 disk structure was presented bv lAndrews et aP ([2009'), 
using the same SMA data. As in that work, we find a moderate cavity size for the DoAr 44 disk, i?cav = 30 AU. 
However, our new model is in much better agreement with the infrared SED, thanks to the inclusion of an inner rim 
and cavity wall in the model structure. To match the relatively faint far-infrared emission, the dust in the outer disk 
is substantially concentrated toward the midplane. Unfortunately, there are no independent constraints on the disk 
inclination or orientation; the adopted inclination could be uncertain by ~20°. The small extent of the 880 /tm emission 
suggests that the disk is quite small (the ring is narrow) , with Rc — 25 AU. While the western peak of the 880 /tm 
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emissi on appears slightly brighter, the model residuals show that the asymmetry is not significant. lEspaillat et alj 
(|2010D modeled the DoAr 44 SED and derived a slightly larger cavity size (36 AU), but within the uncertainties. 

LkCa 15 Located in the Ta urus star-forming r egion, this transition disk has been studied extensively since its cavity 

was first resolved at 1.3 mm by lPietu et al.l (|2006l ). Because of the large spatial extent and prominent depression in the 
880 /im emission, these new SMA observations make for a spectacular example of a resolved transition disk image. The 
cavity size derived here, i?cav = 50 AU, is in excellen t agreement with the initial measurement (46 AU) by IPietu et al.l 
(gOOe). The SED modehng by Esp aillat et all (|2010[ ) argues for a slightly larger cavity (58 AU), although that would 
be difficult to reconcile with the well-defined 880 /im visibility null (see Figure [71). The SEDs for the LkCa 15 and DoAr 
44 disks are quite similar, with both calling for relatively settled vertical dust distributions. An examination of the 
modeling results in this case reveals a striking discrepancy with the data: the model predicts far too little millimeter 
emission inside the cavity (leaving 6 a residuals at the stellar position) . The integrated residual flux density is small, 
^lOmJy, but potentially important because it can provide a more direct estimate for the dust mass inside the cavity. 
For the sake of homogeneity with the rest of this sample, we postpone a more detailed modeling of the LkCa 15 disk 
cavity until higher resolution images are available. 

RX J1615-3255 Perhaps less well-known t han th e ot her d isks in this sample, this system has been kinematically 

tied to the Lupus association, with d » 185 pc (jMa karov 2007). Although it had been identified as a weak- line T Tauri 
star (|Hemzelll976l : lKrautter et al.lll997D. infra red observations by Spitzer have only recently claimed RX J1615-3255 as 
a transition disk candidate ( Merm et al.ll2010l ). The new SMA data presented here represent the first spatially resolved 



images of thermal emission from this disk. A bright, rotating CO J=3— 2 emission structure is also detected, providing 
a preliminary estimate of the viewing geometry. Although much of the outer disk emission is spatially filtered in the 
synthesized 880 /im maps shown in Figures [T] and [71 the visibilities clearly indicate that RX J1615-3255 hosts a rather 
large disk, with i?c = 115 AU (in good agreement with the extent of the CO emission). A particularly low-density 
cavity is present out to i?cav = 30 AU. In fact, we were forced to empty the inner disk inside a radius of 0.5 AU to 
avoid producing a near-infrared emission signature above photospheric levels (therefore, no inner rim was employed 
in this model; see Table HJ. To reproduce the comparatively faint infrared excess beyond the dip in the SED, the 
scale heights were decreased to values consistent with substantial dust settling. The resulting low dust temperatures, 
coupled with the large disk size, lead to a relatively high estimate for the total disk mass (~O.13M0, almost 12% of 
the stellar mass). 

GM Aur This star in the Taurus star- forming region is perhaps the canonical example of a transition disk , with 

many detailed discussions o f the implications of its distinctive infrared SED (e.g. JStrom et al.l[T989: Skrutskie et al.l 
199a iKoerner et all [1991 iChiang fc GoldreichI [l999l: iRice et al.l l2003aD . Using some of the same SMA data. 



Hughes et al.1 ( 200931) first resolved the central dust cavity. We include new 880 ^m SMA observations here, with 



short (C) and intermediate (E) antenna spacings that provide a substantially imp roved Fourier samplin g near the visi- 
bility null. The cavity size derived here, i?cav — 28 AU, is larger than the preferred iHughes et al.l (|2009a[ ) value (20 AU) 
because the new SMA observations provide a more accurate location for the 880 /xm visibility null (at a slightly shorter 
deprojected baseline). However, the alternative disk model discussed bv IHughes et al.l (|2009al ) employed a larger cavity 
(26 AU) that is in good agreement with our model. The long-baseline SMA data is quite iioisy in this case; improved 
sensitivity would help provide a more accurate cavity size estimate. Like the IHughes et al.l (|2009aD results, our model 
successfully matches a high resolution 1.3 mm dataset from the Plateau de Bure interferometer. The GM Aur disk is 
also found to be quite large, with a characteristic size Rc — 120 AU (although most of that emission has been filtered 
out in the images displayed here; see Table [JJ . 

DM Tau Although the gas disk around this cool, young star in the Taurus clouds has been studied extensively 

(e.g., iGuilloteau fc Dutrevi ri998: Dartois et al. 2003; Pictu et al. 2007), our new SMA observations permit the first 
sub-arcsec ond resolutio n exa mination of the 880 /xm dust continuum emission. Based on the shape of its Spitzer IRS 
spectrum, iCalvet et al.l (|2005i) suggested that the DM Tau disk has a relatively small cavity of radius ~3 AU. Since 
such a small cavity cannot be resolved on the longest SMA baselines, it was a surprise to measure a clear 880 /xm 
emission ring and visibility null indicative of a much larger dust-depleted cavity, i?cav = 19 AU. Reconciling this large 
cavity size in the SMA data with the morphology of the IRS spectrum was a challenge, and our current model should 
be considered preliminary until the latter data is better matched. To account for the absence of an infrared excess 
shortward of ^8 /im, we were forced to remove all dust inside a radius of 1 AU (there is no inner rim in this case; see 
Table HI). Like the cases of GM Aur and LkCa 15, we find that the DM Tau disk is rather large, with Rc = 135 AU, 
in good agreement with the extended CO emission studied elsewhere. 

WSB 60 The dust-depleted cavity in the disk a round this very cool star (spectral type M4) in the core of the 

Ophiuchus star- forming region was first discovered bv [Andrews et all (|2009D . Compared to that initial study, we find a 
slightly smaller cavity here (i?cav = 15 AU, compared to 20 AU), partly related to our choice to fix the surface density 
gradient to 7 = 1 in the new models. The cavity is just barely resolved with our SMA observations, and it shows 
no obvious signatures of dust depletion at infrared wavelengths. The latter point drives us to infer the smallest dust 
depletion factor inside the cavity for this sample, corresponding to only a factor of ~50x lower densities in the inner 
disk than the nominal continuous model. Considering the "normal" appearance of the infrared SED, it is natural to 
speculate that there are potentially many other such disks that have similarly large cavities but will only be identified 
when high resolution radio images are available. 
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Table 5 

Sources for SED Data 



Name 


optieal 


infrared 


Spitzer IRS 


(sub)-mrn 


(1) 


(2) 


(3) 


(4) 


(5) 


iVi vv / oo 




12 3 




1 9 


SAO 206462 


2 


1, i, 5 


2 


3, 4, 5 


LkHo 330 


3 


1, 3, 6, 7 


3 


5, 6, 7, 8 


SR 21 


4 


1, 3, 6 


4 
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